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In  the  second  phase  of  the  program,  new  architectures  were  developed 
for  signal  processing  using  time  integration  and  a  comparative  analysis  of 
alternate  techniques  for  one-dimensional  spectrum  analysis  was  performed. 
Coherent  and  non-coherent  time  integrating  schemes  were  considered  and  a  new 
hybrid  space/time  integrating  approach  was  detailed.  The  hybrid  approach 
was  shown  to  exhibit  important  characteristics  related  to  sensitivity  and 
j  dynamic  range,  including  the  suppression  of  effects  due  to  optical  scattering, 
special  sidelobes  and  multisignal  handling  ability.  A  new  approach  to 
ambiguity  function  processing  that  utilizes  space  and  time  integration  in 
two-dimensions  exhibits  superior  performance  in  dynamic  range  and  permits 
high  resolution  Doppler  search. 

The  final  phase  of  the  program  continued  the  extension  of  optical 
architectures  development  that  resulted  in  several  new  space/time  integrating 
system  configurations  for  correlation  and  triple  product  processing  and  a 
space-integrating  triple  product  and  four  product  ambiguity  processor. 
Experimental  results  were  presented  to  demonstrate  the  application  of  chirps 
as  spatially  distributed  local  oscillators  for  interferometric  spectrum 
analysis  and  to  provide  a  quantitative  evaluation  of  a  time-integrating 
ambiguity  function  processor. 
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INTRODUCTION. 


1. 


This  program  was  a  three  year  effort  to  analyze  and  develop 
new  acousto-optic  techniques  for  signal  processing  that  utilize  time 
integration  in  one  and  two  dimensions.  The  following  paragraphs  in 
this  section  summarize  the  results  of  the  first  two  years  (phases  1  and 
2) .  Sections  2  and  3  present  a  review  of  the  results  of  the  final 
phase,  which  includes  the  extension  of  system  architectures  to  achieve 
high  performance  processing  and  an  investigation  of  realizability 
considerations,  including  laboratory  measurements  of  some  time-integrating 
processing  systems. 

1. 1 _ Summary  of  Phase  I . 


1.1.1  Architectures. 


Time  integrating  optical  signal  processing  has  been 
generalized  for  one-  and  two-dimensional  complex  signal  processing. 
Signal  processing  architectures  utilizing  actively  generated  reference 
waveforms  such  as  chirps  realize  a  wide  range  of  algorithms,  and 
variable  time  integration  allows  further  flexibility.  The  class  of 
signal  processing  operations  achieved  with  time  integrating  acousto¬ 
optic  techniques  has  been  described,  and  hybrid  architectures  using 
both  time  and  space  integration  have  been  discussed.  Particularly 
attractive  is  the  multi-purpose  capability  of  the  time  integrating 
optical  processor,  e.g.,  the  same  optical  system  is  used  for  both 
spectral  analysis  and  ambiguity  function  processing. 

1.1.2 _ Optical  Processor  Realizations. 

Several  optical  realizations  were  analyzed  including  inter¬ 
ferometric  and  non-interferometric  implementations.  Interfereometric 
techniques  are  based  on  linear  electric  field  modulation  and  non- 
interferometric  techniques  are  based  on  linear  intensity  modulation. 
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1.1.2 


The  concept  of  interferometric  detection  with  electronic  reference  was 
introduced.  In  this  implementation,  a  local  oscillator  is  added  to  the 
acoustic  modulation  rather  than  a  coherent  optical  reference  added  to 
the  image.  This  method,  therefore,  permits  non-coherent  optical 
implementation,  with  the  advantage  of  directly  modulated  diode  light 
sources  and  increased  immunity  to  artifacts  of  coherent  optical  imaging 
systems.  The  correlation  is  performed  at  a  carrier,  thereby  enabling 
complex  operation.  Furthermore,  the  interferometric  method  circumvents 
the  difficult  requirements  for  acoustic  modulation  at  a  high  diffraction 
efficiency  bias  point,  which  is  necessary  for  linear  operation  in  the 
non-interferometric  technique. 

1.1.3 _ Statistical  Analysis. 

Fundamental  limitations  were  established  with  regards  to 
signal-to-noise  and  signal-to-bias  ratios.  The  reduction  in  signal- 
to-noise  ratio  due  to  detection  noise  and  bias  variation  was  analyzed. 

It  was  shown  that  the  average  bias  reduces  the  dynamic  range  and  that 
the  bias  variation  cause  a  loss  in  sensitivity.  Expressions  were 
derived  for  the  optical  processor  output  that  include  the  signal  depen¬ 
dent  bias  contribution.  This  contribution  may  be  filtered  or  subtracted. 
Several  schemes  were  analyzed. 

The  optical  implementation,  utilizing  acousto-optic  input 
and  integrating  image  sensor  output  devices,  creates  a  processor  that 
is  highly  compatible  with  signal  processing  systems.  Such  implementa¬ 
tion  affords  very  large  time-bandwidth  product  signal  processing  without 
the  input  signal  storage  requirement  associated  with  spatial  integrating 
methods.  The  requirement  for  high  resolution  large  dynamic  range  output 
image  sensors  becomes  the  key  device  limitation.  Acousto-optic  devices 
are  available  with  time-bandwidth  product  much  greater  than  the  number 
of  resolvable  image  samples. 
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The  key  attributes  of  time  integrating  techniques  are 
summarized  in  the  following:  extremely  large  time-bandwidth 
correlations  may  be  performed,  independent  of  the  device  time-bandwidth 
product;  a  flexible,  multi-purpose  processor  is  realized  by  use  of 
actively  generated  reference  waveforms  and  variable  time  integration; 
an  important  class  of  two-dimensional  algorithms,  including  complex 
spectral  analysis  and  ambiguity  function  processing,  may  be  performed 
without  having  to  store  the  entire  time  history  as  a  spatial  history; 
the  system  may  be  implemented  with  non-coherent  diode  sources,  acousto¬ 
optic  devices,  and  integrating  image  sensors. 

1.2 _ Summary  of  Phase  II. 

Efforts  during  this  phase  concentrated  on  system 
architectures  for  spectrum  analysis,  including  a  comparative  analysis 
of  various  approaches  to  one-dimensional  processing,  dynamic  range 
studies  and  implementation  requirements,  and  coherent  hybrid  approaches 
to  ambiguity  function  processing. 

1.2.1 _ One-Dimensional  Spectrum  Analysis. 

Many  acousto-optic  techniques  are  available  for  spectrum 
analysis.  The  interferometric  schemes  can  be  divided  into  three  groups 
time-integrating,  space-integrating  and  hybrid  space/time  integrating. 
These  have  been  analyzed  to  establish  their  capabilities  with  respect 
to  each  other  and  in  comparison  with  an  acousto-optic  (non-interfero- 
metric)  power  spectrum  analyzer  (which  we  will  call  a  Bragg  receiver) . 

A  Bragg  receiver  measures  the  signal  power  spectrum,  so  that 
long  detector  integration  can  be  used  to  gain  signal-to-noise  ratio  for 
narrowband  signals  and  to  ease  the  detector  read  rate  requirements.  A 
major  difficulty  with  the  Bragg  receiver  is  the  limitation  on  the 
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dynamic  range  imposed  by  the  detector  noise.  A  common  characteristic 
of  all  the  interferometric  approaches  is  the  considerable  improvement 
of  the  dynamic  range  attained  by  measuring  the  complex  amplitude 
spectrum. 

1.2. 1.1  Time  Integrating  Techniques. 

A  number  of  configurations  for  time  integration  spectrum 
analysis,  both  coherent  and  incoherent,  have  been  analyzed.  Their  major 
characteristics  may  be  summarized  as  follows: 

•  The  bandwidth  is  limited  to  about  25%  of  the  Bragg 
cell  bandwidth  by  chirp  transition  effects. 

•  Chirps  that  are  not  highly  uniform  in  amplitude  and 
linear  in  phase  can  produce  spurious  signals.  A 
noisy  chirp  will  produce  a  noisy  spectrum  when  a  signal 
is  present. 

•  Multiple  signal  handling  capacity  is  low.  Unit 
modulation  can  occur  only  for  a  narrowband  signal. 

The  dynamic  range  therefore  degrades  as  the 
spectrum  becomes  more  dense. 

•  Counterpropagating  schemes  sometimes  require  an 
external  delay  of  t  (equal  to  the  Bragg  cell  delay) 
for  the  chirp  signal  applied  to  one  of  the  cells. 

Also,  whenever  a  baseband  signal  is  mixed  with  a  CW 
or  chirp,  the  product  must  be  generated  single  side¬ 
band  to  avoid  a  3  dB  loss  in  dynamic  range. 

•  Optical  scatter  noise  is  suppressed.  Light  scattered 
from  optical  elements  is  not  Doppler-shifted,  and 
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1.2. 1.1 


therefore  the  only  components  in  the  exposure  that  do 
not  integrate  to  zero  are  products  containing  only 
scatter  elements.  The  dynamic  range  referred  to 
scatter  noise  is  twice  as  large  (in  dB)  for  time- 
integrating  techniques  as  it  is  for  a  Bragg  receiver. 

1.2. 1.2  Space  Integrating  Interferometric  Techniques. 

The  conf iguration  considered  here  is  one  using  a  pulsed 
coherent  source,  such  as  a  laser  diode,  to  generate  a  snapshot  of  the 
amplitude  spectrum.  The  chief  attraction  of  this  technique  is  the 
removal  of  the  requirement  for  a  chirp.  An  additional  favorable  aspect 
of  the  pulsed  laser  scheme  is  that,  unlike  the  time-integration 
techniques  and  wider-band  hybrid  techniques,  the  reference  wave  does 
not  need  to  be  diffracted  by  a  Bragg  cell.  Since  the  reference  wave 
needs  to  illuminate  all  the  diodes  in  the  detector  array,  a  low 
efficiency  diffraction  of  the  reference  wave  could  have  a  serious 
impact  on  the  system  light  budget,  especially  at  wide  bandwidths . 

Optical  scatter,  however,  is  not  suppressed  unless  a  supple¬ 
mentary  Bragg  cell  is  used  in  the  reference  arm  of  the  interferometer, 
complicating  the  design. 

Laser  diodes  are  presently  the  only  realistic  pulsed  light 
source  for  this  application.  Diode  sources  with  adequate  coherence 
and  sufficient  peak  and  average  output  powers  are  available,  but  the 
transverse  mode  structures  of  these  devices  are  grossly  inadequate. 

1.2. 1.3  Coherent  Hybrid  Techniques. 

The  coherent  hybrid  approach  to  interferometric  spectrum 
analysis  described  here  uses  both  space  integration  and  time  integration. 
In  summary 
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•  Spectral  sidelobes  are  down  by  a  factor  of  two  in  dB 
compared  with  the  above  techniques  and  with  the  Bragg 
receiver.  Sidelobes  of  the  very  bright  zero  order 
(undiffracted  by  the  acoustic  wave)  that  extend  into 
the  signal  band  are  similarly  suppressed.  This 
inherent  sidelobe  suppression  is  the  result  of  the 
combined  effect  of  the  localization  of  the  signal 
illumination  at  the  detector  by  spatial  integration 
(i.e.,  the  optical  Fourier  transform)  and  localization 
of  the  interference  fringe  pattern  by  time  integration. 

•  As  in  the  purely  time-integration  techniques,  scatter 
noise  from  optical  elements  is  suppressed.  Dynamic 
range  relative  to  scatter  noise  is  doubled  in  dB 
compared  with  that  of  the  Bragg  receiver. 

•  Hybrid  techniques,  unlike  time-integrating  techniques 
are  not  sensitive  to  chirp  noise  and  chirp-related 
intermodulation  distortion. 

•  The  multiple  signal  handling  capability  of  a  hybrid 
spectrum  analyzer  is  not  degraded  from  that  of  a  Bragg 
receiver.  Whereas  time-integration  techniques  degrade 
in  dynamic  range  by  10  log  N  due  to  the  buildup  of  bias 
terms  when  N  equal  amplitude  signals  are  present, 
hybrid  and  space-integrating  techniques  spatially 
separate  the  light  associated  with  signals  of  different 
frequency. 

1.2.2 _ Dynamic  Range  and  Implementation  Factors. 

Dynamic  range  for  interferometric  detection  was  shown  to  be 

,  2 
limited  by  the  detector  noise  (of  variance  a  )  to 
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DR  =  (Vgat/o)2/8  , 

where  Vgat  is  the  detector  saturation  voltage.  Since  for  a  Bragg 
receiver  DR  =  V  t/o,  the  gain  in  dynamic  range  is  considerable. 

Shot  noise  due  to  the  bias  can  be  significant  in  large 
dynamic  range  systems.  Large  saturation  values  for  the  detector  and 
high  light  powers  are  required  to  avoid  this  limitation. 

Other  serious  dynamic  range  limits,  as  mentioned  above,  can 
be  overcome  with  the  correct  system  architecture.  Time-integration 
suppresses  scatter  noise,  and,  in  addition,  hybrid  techniques  overcome 
chirp  distortion  and  noise  and  also  suppress  sidelobes. 

Chirp  linearity  requirements  and  the  effects  of  nonlinearities 
were  investigated.  It  was  found  that  resolution  considerations  permit 
larger  percent  linearity  error  for  time  integrating  techniques,  but  the 
requirement  to  suppress  intermodulation  terms  imposes  a  much  more  severe 
linearity  requirement  for  time  integrating  than  for  hybrid  techniques. 

1.2.3 _ Ambiguity  Function  Processing. 

A  hybrid  ambiguity  function  concept  was  developed  that 
utilizes  a  space-integrating  ambiguity  function  generator  in  one  arm  of 
an  interferometer  and  a  chirp-generated  spatially  distributed  local 
oscillator  in  the  other  arm.  Time  integration  detection  yields  an 
estimate  of  the  ambiguity  function  with  finer  Doppler  resoltuion  (than 
a  space  integrating  technique)  but  with  a  signal-dependent  bias 
localized  in  the  delay  direction.  It  was  shown  that  a  multiplicity  of 
correlation  peaks  does  not  result  in  a  decreased  dynamic  range  due  to 
bias  buildup,  as  it  does  in  purely  time  integrating  techniques,  provided 
that  the  peaks  are  separated  by  at  least  one  delay  resolution  element. 


1-7 


ARCHITECTURES. 


*s 


2. 


It  has  become  apparent  that  there  is  an  innumerable  variety 
of  architectures,  implementation  variations  and  algorithms  possible 
with  acousto-optic  techniques ,  and  that  optimum  architectures  for  one 
algorithm  are  not  generally  optimum  for  another  algorithm.  Further¬ 
more,  ithplementational  variations  may  differ  widely  in  performance, 
relative  to  signal-to-noise  ratio,  bandwidth  and  other  important 
parameters.  It  is  further  recognized  that  the  driving  motivation  for 
A-0  processing  is  the  ability  to  process  data  at  very  high  rates.  With 
these  considerations  in  mind  it  is  highly  relevant  to  compare  and 
contrast  A-0  techniques,  time  integrating,  space  integrating  and  hybrid 
architectures  for  correlation,  generalized  triple  product  processing 
and  ambiguity  function  generation.  In  the  paragraphs  below  we  discuss 
some  extensions  of  previously  described  configurations  and  introduce 
several  new  architectures. 

2.1. _ Correlation. 

In  this  section  one  dimensional  A-0  techniques  for  signal 
correlation  are  discussed.  The  performance  of  time-integrating, 
space-integrating  and  hybrid  configurations  are  compared  with  one 
another  and  with  respect  to  the  performance  of  an  ideal  correlator. 

2.1.1  Ideal  Correlator. 


Let 

S^t)  =  S(t-tQ)  +  nx(t)  (2-1) 

S2(t)  =  S(t)  +  n2(t) 

be  complex  signals  plus  noise;  the  cross-correlation  function  between 
these  signals  is  defined  as 
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where  R^(0),  1^(0)  and  Rg(0)  are  the  input  noise  and  signal  powers, 
respectively. 


1 


2.1.2 _ Acousto-Optic  Time-Integrating  Correlator. 

There  are  many  variations  of  A-0  time-integrating  correlators. 
For  the  present  discussion  we  will  choose  the  noncoherent  configuration 
in  which  a  light  source  is  modulated  in  intensity  according  to 


I  ^  ( t )  =  1  +  m^  |  (t)  |  cos  [27rf  t  +  <p^(t)]  , 


(2-7) 


and  the  Bragg  cell  is  driven  so  that  the  diffraction  efficiency  at  a 
point  x  from  the  transducer  is  I2(t-r),  where 


I2(t)  =  1  +  m2  |  S2  (t)  |  cos  [2irft  +  4*^  (t)  ] 


(2-8) 


It  is  assumed  that  f  3B/2.  The  input  gains  and  m2  are  adjusted 
so  that 


m2  Etls^t)  |2] 


m2  Et 1 S2 (t)  |2] 


«  m^[Rs(0)  +  Rx(0)]  <  1 


=  m2[Rs(0)  +  R2(0)]  <  1 


(2-9) 


The  light  diffracted  from  the  Bragg  cell,  imaged  onto  the  detector 
array,  is  integrated  for  timeT,  after  which  the  voltage  output  from 
the  diode  at  t  is 


V(t)  = 


f  H 

j  C\ 


(t)I-(t-T)dt  +  V,(T) 


(2-10) 


where  V_  is  the  accumulated  bias  and  V-  is  the  electrical  noise  added 
B  d 

by  the  detector.  From  Equations  2-7  and  2-8, 


V  (t )  =  VB  +  Vg(T)  +  Vd(T) 


(2-11) 
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where 


VS(T) 


and  where 


mlm2VB 

- 2 - 


|  2^)1  008  f 2irf T 


+  (T  )  ] 

i 


R1  2(x)  *  ^Rl,2  I  exp  I  j  (r )  ] 


(2-12) 


With  appropriate  demodulation  the  signal  can  be  detected  with  an  output 
signal-to-noise  ratio  given  by  L 


where  is  the  mean  of  | Vd (t ) |  times  the  number  of  detector  elements 
per  delay  bin.  With  the  signals  of  Equations  2-1  and  2-2,  this  reduces 
to 

SNR(tQ) 

SNR' (tQ)  =  1  +  aBT  .  (2-14) 

M2 

In  this  example  a  =  8.  SNR(tg)  is  the  SNR  for  an  ideal  correlator 
(Equation  2-5),  and  M  =  Vg/crd  represents  the  dynamic  range  of  the 
detector.  A  value  of  equal  to  the  saturation  voltage  determines 
the  limiting  performance.  Equation  2-14  quantifies  the  limit  of 
processing  gain  for  a  noncoherent  time-integrating  correlator.  For 
example,  if  the  saturation-to-rms  noise  of  an  element  of  a  detector 
array  is  1000,  and  if  there  are  8  sampling  elements  per  correlation 
width,  the  limiting  value  of  is  8  x  10  ,  and  60  dB  of  processing 
gain  is  achievable.  This  may  not  be  adequate  for  some  applications 
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when  used  as  an  auto-correlator  with  a  low  input  SNR.  The  dynamic 
range  limitation  is  aggravated  by  the  presence  of  other  signals  in  the 
signal  band,  which  contribute  to  the  bias. 

The  bandwidth  of  the  Bragg  cell  in  this  noncoherent  time- 
integrating  configuration  must  be  at  least  twice  the  signal  bandwidth, 
since  a  reference  tone  must  be  placed  at  a  frequency  at  least  3B/2 
away  from  the  center  of  the  band. 


Noncoherent  schemes  using  two  Bragg  cells  permit  the  use  of 
a  bright  light  source,  such  as  an  arc  lamp,  that  cannot  be  directly 
modulated.  Two  Bragg  cells,  one  imaged  onto  the  other  in  a  counter- 
propagating  mode,  also  permit  correlations  to  be  measured  over  twice 
the  relative  delay  as  a  single  cell  correlator.  The  same  bandwidth 
and  dynamic  range  considerations  as  outlined  above  apply,  however. 


A  coherent  time-integrating  correlator  can  be  implemented 
so  that  the  Bragg  cells  have  the  same  bandwidth  as  the  signal.  We 
assume  a  configuration  in  which  two  Bragg  cells  are  in  separate  arms 
of  an  interferometer,  such  that  the  signal  amplitudes,  represented 
by  the  optical  amplitudes  diffracted  from  the  cells,  are  imaged  onto 
a  linear  detector  array.  The  power  density  on  the  array  is 

I  ( t , x )  =  im1S1(t-x)  +  m2S2  (t+T)ej2lrfT|2  ,  (2-15) 


where  the  exponential  describes  the  effect  of  the  offset  angle 
between  the  beams .  The  time  integrated  array  output  is 


V  (t  ) 


2m^m  2 


|  R12  (t)  |  cos  [2TTfT  +  4>  (t )  ] 


+  V  - 


(2-16) 
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where  the  bias  is 


VB  =  [  |S1(t-x) |2dt  +  ^  f  |s2(t+T)|2dt 

J  n  J  r\ 


(2-17) 


It  can  be  shown  that  the  output  SNR  for  this  coherent  technique  used 
as  an  auto-correlator  is  given  by  Equation  2-14  with  a  =  2,  indicating 
that  the  processing  gain  limitation  is  6  dB  less  serious  than  for  a 
noncoherent  technique.  Interference  can  degrade  the  performance  of  a 
coherent  correlator  also. 


2.1.3 _ Noncoherent  Space-Integrating  Correlator. 

The  noncoherent  correlator  with  two  Bragg  cells  can  be 
used,  with  slight  modification,  as  a  space-integrating  processor.  If 
the  input  to  one  cell  is  a  time-reversed  version  of  the  input  to  the 
other,  and  if  the  linear  detector  array  is  replaced  by  a  single  large- 
area  wideband  detector,  the  cross-correlation  of  the  signals  is 
contained  in  the  time-varying  output  of  the  detector: 

^1^0  /v 

V(t)  =  VB  +  -  j-  -  R12(2t)cos[4irft  +  ip  (2t)  ]  +  Vd  (2-18) 

The  correlation  function  is  on  a  temporal  carrier  and,  therefore, 
requires  a  detector  with  a  wide  bandwidth.  A  detector  with  a  large 
active  area  has  a  large  capacitance,  making  it  difficult  to  provide 
wideband  response  with  high  dynamic  range.  A  space-integrating 
correlator  provides  processing  gain  equal  to  the  time-bandwidth  product 
of  the  Bragg  cells,  and  therefore  limits  its  usefulness  when  large 
processing  gain  is  required,  particularly  in  an  auto-correlator.  When 
used  as  a  cross-correlator,  the  frequency  of  both  signals  must  be  equal 
to  an  accuracy  of  the  reciprocal  of  the  integration  time  (as  for  all 
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linear  correlators) .  A  sequential  search  through  frequency  may  be 
required;  to  minimize  total  search  time  it  is  desirable  to  provide  the 
minimum  processing  gain  required  to  detect  the  signal.  In  those 
cases  where  space  integration  provides  adequate  gain,  a  space-integrating 
correlator  may  be  a  useful  device  for  rapid  search.  Coherent 
implementations,  however,  can  provide  significant  advantages.  These 
advantages  include  the  capability  for  parallel  search  over  frequency; 
such  configurations  perform  simultaneous  cross-correlation  and  spectrum 
analysis  and  are  referred  to  as  ambiguity  function  processors. 


2.1.4 


Coherent  Space/Time-Integrating  Correlator. 


The  technique  described  here  is  a  coherent  optical  approach 
that  exploits  the  advantages  of  space  integration  processing  and  yet 
permits  the  extended  gain  provided  by  time-integration.  It  is  also 
hybrid  in  the  sense  that  the  technique  includes  a  digital  processing 
step;  this  occurs,  however,  after  time  integration  where  the  data  rates 
are  compatible  with  digital  techniques. 

The  approach  is  illustrated  in  Figure  2-1,  and  is  basically 
an  integrating  acousto-optic  power  spectrum  analyzer  followed  by  a 
digital  Fourier  transform.  The  input  to  the  Bragg  cell  is  the  sum  of 
the  two  signals  being  correlated.  The  optical  power  integrated  over 
time  T2  on  the  detector  results  in  an  exposure 


E  { f )  = 


A1(f,t)  +  A2(f,t) 


(2-19) 


where  for  A^  and  A2  are  the  optical  Fourier  transform  estimates  of 
signals  and  S2 # 
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Figure  2-1.  Space/Time  Integrating  Correlator. 
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A  ( f ,  t ) 


S (t-i ) e 


j27rf 


dt 


(2-20) 


where  is  the  Bragg  cell  delay.  The  detector  output  is  proportional 
to  E(f),  plus  detector  noise  V^fx).  The  Fourier  transform  of  the 
result,  performed  digitally,  is 


W  (x ) 


(x  )  |  cos  'j  (x) 


+  n(x) 


(2-21) 


The  FT  of  the  bias  term  occurs  at  x  =  0  and  has  been  ignored.  Vg  is 
the  bias  level  of  the  detector,  n(x)  is  the  FT  of  the  detector  noise, 
where 

E  [n2  (x  )  ]  =  (B/T1)  E  [V^j  (f )  ] 

and 

P 

=  [Rg^O)  +  Rg2(0)]/B  .  (2-22) 


Equation  2-21  indicates  that  the  real  output  of  the  correlation 
estimate  is  obtained.  A  second  measurement  with  S2  shifted  in  phase  by 
90°  yields  the  quadrature  component. 

A  calculation  of  the  output  signal-to-noise  ratio  for  signals 
and  noise  are  defined  in  Equation  2-1  leads  to 
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SNR’ (tQ) 


1  + 


SNR(tg) 


2 


(2-23) 


i.e.,  a  =  2/BT^ .  The  processing  gain  achievable  for  a  hybrid  space/ 
time  integrating  correlator  is,  therefore,  larger  than  that  for  a 
purely  time-integrating  correlator  by  at  least  a  factor  BT^.  This 
improvement  in  performance  over  purely  time-integrating  techniques  is 
directly  analogous  to  the  signal  handling  ability  and  dynamic  range 
improvement  in  the  hybrid  space/time  spectrum  analyzer  and  the  hybrid 
ambiguity  function  processor  reported  earlier  (see  the  Interim  Technical 
Report,  1980) .  This  technique  of  measuring  the  joint  spectral  density 
with  the  detector  array  also  permits  the  excision  (or  at  least 
limitation)  of  narrowband  interference  signals  occurring  in  the  band 
before  digital  Fourier  transformation.  It  is  also  tolerant  of  point 
blemishes. 


2 . 2 _ Triple  Product  Correlator. 


A  triple  product  processor  performs  the  estimate 


R(x^, t2  ) 


S1(t)S2(t-x1)S3(t-x3)dt 


(2-24) 


which  can  be  implemented  optically  in  a  variety  of  ways.  Below  we 
contrast  two  techniques,  one  purely  time  integrating,  the  other  a 
scheme  that  employs  both  space-integration  and  time-integration. 

2.2.1 _ Noncoherent  Time  Integrating  Processor. 

The  noncoherent  correlator  of  subsection  2.1.2  can  be 
expanded  to  two  dimensions  by  passing  the  light  diffracted  from  the 
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Bragg  cell  through  a  second  cell  orthogonal  to  the  first  driven  by 
the  third  signal.  The  second  Bragg  cell  produces  a  diffracted 
intensity. 


I3 (t-x  2>  =  1  +  m3s3(t-T3) 


(2-25) 


S3  is  assumed  to  be  real.  Both  Bragg  cells  are  imaged  onto  a  two 
dimensional  photodetector  array.  The  important  term  of  the  triple 
product  exposure  integral  over  time  T  is 


VS(t1't2)  = 


m  1  m?m ^Vp 


A 

[  R(t1#t2)  I  cos  [2irf t  1+'p  ( T 1 » T 2 )  ]  (2-26) 


where  the  correlation  estimate  is  on  a  carrier  in  the  x^  direction,  in 
this  example.  The  output  SNR  can  be  shown  again  to  be  given  by 
Equation  2-14  with  <*  =  8.  This  means  that  the  noncoherent  two- 
dimensional  triple  product  correlator  permits  the  same  processing  gain 
as  the  one-dimensional  device.  For  cross-correlation  of  noisy  signals, 
where  large  processing  gain  is  desired,  this  technique  will  frequently 
be  inadequate,  especially  if  interfering  signals  are  present.  The 
Bragg  cell  bandwidth  again  must  be  at  least  twice  the  signal  bandwidth. 
(The  full  Bragg  cell  bandwidth  may  be  used  if  the  measurement  is  made 
without  a  carrier  twice.)  Coherent  time-integrating  techniques  can 
provide  some  improvement,  as  in  the  one-dimensional  case,  but  not  nearly 
as  much  as  hybrid  space/time  integrating  techniques. 


2.2.2 


Coherent  Space/Time  Processor. 


The  optical  configuration  shown  in  Figure  2-2  generates  a 
two  dimensional  cross-spectral  density  that  is  transformed  into  the 
triple  product  integral  by  a  digital  Fourier  transform.  Light  from  a 
coherent  source  is  split  into  two  arms  each  containing  a  Bragg  cell. 
The  light  diffracted  from  each  of  the  orthogonally  oriented  cells. 
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driven  by  S2  and  is  optically  Fourier  transformed  onto  a  two- 
dimensional  detector  array,  where  the  optical  power  density  is 


I(f1,f2,t) 


|A1(f1,t)  +  A2(f2,t) 


(2-27) 


where  and  A2  are  the  spectral  amplitude  estimates  of  and  S2  as 
defined  in  Equation  2-20.  The  light  source  is  modulated  in  intensity 
by  a  third  Bragg  cell.  The  time  integrated  power  density  produces  a 
detector  output  signal 


v(fi,f2) 


V. 


B 


[1  +  S3  (t)  ]  |  Ax  (f  1#t) 


(2-28) 


+  A2(f2,t)l2dt 


The  two  dimensional  digital  Fourier  Transform  produces 
auto-correlation  lines  at  =  0,  t2  =  0,  a  diagonal  cross-correlation 
line  between  and  S2,  and  the  term  of  interest 

Vg(t)  =  VB| R(t1,t2) | cos^ (t1,t2)  (2-29) 

The  second  quadrature  component  is  obtained  by  a  second  measurement 
with  the  identical  input  signals  but  with  one  of  the  signals  having  a 
90°  phase  shift  imposed.  The  processing  gain  permitted  with  this 
technique  is  limited  by  Equation  2-14  with  a  =  2/BT^,  as  in  the  one¬ 
dimensional  space/time  integrating  correlator.  In  addition  to  the 
large  gain  that  is  achievable,  this  technique  permits  relatively  simple 
suppression  of  narrowband  interference  signals  after  time-integration 
but  before  the  digital  transform.  Point  blemishes  on  large  two- 
dimensional  photodetector  arrays  are  almost  impossible  to  avoid  with 
the  present  state  of  technology,  and  these  are  particularly  troublesome 
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when  they  occur  in  the  correlation  domain.  However,  in  the  spectral 
domain,  as  in  the  present  case,  blemishes  are  not  particularly  damaging, 
as  their  effect  is  more  or  less  uniformly  distributed  over  the 
correlation  field  after  Fourier  transformation. 

2.2.3 _ Space-Integrating  Triple  Product  Processor. 

A  space-integrating  processor  that  computes  the  triple 
product  integral  is  shown  in  Figure  2-3.  The  first  Bragg  cell  has  a 
diffracted  optical  amplitude  that,  when  imaged  onto  the  second  Bragg 
cell,  is  proportional  to  the  signal  traveling  opposite  to  that  in  the 
second  Bragg  cell.  The  light  diffracted  from  the  second  cell  is 
(x+t) S2 (x-t) .  This  is  imaged  in  the  x  direction  through  the  third 
Bragg  cell  onto  plane  P;  the  plane  of  the  third  cell,  driven  by  S3(t), 
is  imaged  in  the  y  direction  onto  plane  P.  The  result  is  Fourier 
transformed  in  the  direction  v  =  (x+y)//2f  and  imaged  in  direction 
2  =  (x-y)//7  onto  the  detector  plane.  The  resulting  amplitude  is 

E(t,w)  =  /  S1[(v+w)//2+t]S2[(v+w)//Z-t]S3[(v-w)//?+t]ej2lTVXdv 

=  ej2lT(w_/7t)zys1(5)S2(C-2t)S3(C-w/I)ej2lT/25zd  .  (2-30) 

At  z  =  0  the  detected  intensity  is 

I E | 2  =  | R[2t,/2w] | 2  ,  (2-31) 

recognized  as  the  square  of  the  triple  product  integral,  with  =  2t 
T 2  =  ^w-  A  linear  photodiode  array  with  parallel  outputs  placed  along 
z  =  0  permits  very  high  speed  processing  in  those  cases  in  which  a 
processing  gain  equal  to  the  time  bandwidth  product  of  the  Bragg  cell 
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is  sufficient.  Detector  processing  considerations  are  similar  to 
those  of  the  one-dimensional  space-integrating  ambiguity  processor. 

A  triple  product  processor  can  be  used  in  a  linear  carrier- 
dependent  mode,  where  frequency  must  be  known.  The  processor  of 
Figure  2-3  can,  in  principal,  be  expanded  to  a  three-product  ambiguity 
processor  by  implementing  a  two-dimensional  detector  array  in  which  all 
elements  have  parallel  outputs.  Very  high  processing  speeds  would  be 
achievable  with  this  processor.  The  detector  array  problem,  however, 
is  extremely  difficult.  Some  compromises  might  be  made  to  provide 
some,  but  not  all,  of  the  frequency  coverage,  such  as  placing  a  number 
(say  N)  of  linear  arrays  equally  spaced  in  the  plane,  with  a  consequent 
reduction  in  Doppler  search  time  by  a  factor  N. 

2 . 3 _ Ambiguity  Function  Processing. 

When  a  correlation  is  performed,  for  example,  between  an 
unknown  signal  and  a  hypothesized  signal,  the  frequency  of  the  unknown 
signal  must  be  known  to  an  accuracy  of  at  least  the  reciprocal  of  the 
correlation  integration  time.  If  it  is  not,  a  search  through  frequency 
is  required.  Cross-correlation  between  a  received  signal  and  a 
multipath  component  also  requires  a  search  through  frequency  if  the 
relative  Doppler  shift  is  not  known  to  the  required  accuracy.  Ambiguity 
function  processors  permit  simultaneous  correlations  to  be  performed 
between  a  signal  of  unknown  frequency  and  a  signal  with  many  hypothe¬ 
sized  frequencies. 

Below  we  discuss  several  types  of  ambiguity  functions  that 
promise  to  be  of  value  in  substantially  reducing  processing  time  in 
applications  requiring  search  through  delay  and  Doppler. 
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2.3.1  _ Triple  Product  Ambiguity  Processor. 

Any  processor  that  generates  the  triple  product  integral 
can  be  used  for  ambiguity  processing  by  making  one  input  a  chirp  and 
premultiplying  one  of  the  signal  inputs  with  the  same  chirp.  For 
reasons  discussed  above  and  in  earlier  reports ,  time-integrating 
processors  have  some  disadvantages  compared  with  space/time  integrating 
hybrid  techniques.  A  difficulty  with  a  wideband  ambiguity  processor 
using  time  integration  is  that  a  wideband  highly  linear  chirp  is 
required,  if  a  triple  product  processor  is  used  directly,  or  else  the 
processor  must  be  reconfigured  to  use  a  narrower  bandwidth  Bragg  cell 
with  which  one  could  use  a  digitally  generated  chirp. 

Wideband  analog  chirps  can  be  generated  with  a  VCO,  for 
example,  but  the  linearity  of  such  a  device,  even  with  the  most 
sophisticated  linearization  techniques  presently  available,  is  not 
nearly  adequate.  Highly  linear  surface  acoustic  wave  (SAW)  devices 
with  wide  bandwidths  are  becoming  available  (see  Section  3) .  These 
suffer  the  severe  disadvantages,  however,  of  having  a  fixed  delay 
and  bandwidth  and  cannot  be  fabricated  with  the  very  large  chirp 
time-bandwidth  product  required  for  high  Doppler  resolution. 

2.3.2  _ One-Dimensional  Space  Integrating  Processor. 

An  extension  of  the  one-dimensional  space-integrating 
correlator  can  be  configured  to  permit  parallel  frequency  search. 

Figure  2-4  illustrates  a  correlator  that  performs  parallel  frequency 
testing.  Coherent  light  diffracted  from  a  Bragg  cell  driven  by  S^(t) 
is  imaged  onto  a  second  cell  driven  by  S2(-t)  so  that  the  signals  are 
counterpropagating .  The  light  diffracted  from  the  second  cell  is 
optically  Fourier  transformed  onto  an  array  of  photodiodes  with 
parallel  outputs.  The  optical  power  density  in  the  array  plane  is 


2-17 


2.3.2 


Continued 


I(t) 


S,  (t-x)s!!(-t-T)ej2irfdt 


=  (x12(f,t)|2 


(2-32) 


The  detector  array  senses  the  absolute  sqare  of  the  ambiguity  function, 
where  time  is  the  correlation  axis  and  position  in  the  diode  array  is 
the  frequency  axis.  Since  there  is  no  time-integration  processing, 
this  technique  is  only  useful  when  the  processing  gain  equal  to  the 
time  bandwidth  product  of  the  Bragg  cells  is  sufficient.  In  such  cases 
this  device  could  be  a  valuable  tool  for  very  rapid  search,  since 
ambiguity  surfaces  can  be  generated  at  very  high  rates  due  to  the 
channelization  in  frequency.  Note  that  for  a  signal  bandwidth 
B  =  200  MHz,  for  example,  with  a  signal  length  T  *  20  ps,  a  Bragg  cell 
bandwidth  of  400  MHz  and  time-bandwidth  product  of  500,  a  search  of 
105  ambiguity  surfaces  per  second  results,  each  surface  containing 
8  x  105  delay-frequency  elements.  The  surfaces  would  be  detected  by  a 
linear  array  of  250  photodiodes,  each  with  a  200  MHz  video  bandwidth. 


For  applications  requiring  only  threshold  detection,  post¬ 
detection  processing  architectures  can  be  several  levels  of  complexity. 
The  technique  providing  the  highest  degree  of  sensitivity  is  to  implement  a 
synchronous  video  integrator  on  each  channel,  thus  noncoherently  adding 
successive  ambiguity  planes.  This  post-detection  integration  generates 
further  processing  gain  over  and  above  that  provided  by  the  coherent 
space  integration.  The  size,  speed,  versatility  and  quantity  of  these 
channel  processors  make  this  approach  presently  impractical.  A  realistic 
approach  is  to  threshold  each  channel  and  process  threshold  crossings 
to  detect  signals  while  minimizing  errors.  The  probability  of  false 
alarm  (Pfa)  is  a  function  of  the  threshold  level,  which  affects  sensi¬ 
tivity.  In  the  above  example  a  threshold-to-noise  ratio  of  12  dB 
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(corresponding  to  Pfa  =  2.5  x  10  for  each  delay  bin)  provides  a  Pfa 
of  10-3  over  a  4000  bin  delay  range.  If  each  frequency  channel  is 
revisited  for  confirmation  the  net  P^a  is  10  ^  for  each  frequency 
channel  over  20  ys,  or  one  confirmed  error  every  20  seconds.  The  same 
false  alarm  rate  can  be  achieved  with  a  5  dB  threshold-to-noise  ratio 
if  each  delay  bin  of  each  frequency  channel  can  be  tested  for  confir¬ 
mation.  The  detailed  behavior  and  requirements  for  various  post-process 
schemes  (including  channel  integration  of  threshold  crossings)  need  to 
be  determined. 

In  an  alternate  arrangement,  the  second  Bragg  cell,  driven 
by  a  digitally  stored  and  time-reversed  signal,  is  narrower  bandwidth 
but  with  larger  delay  than  the  first  cell,  driven  by  a  wideband 
signal.  The  imaging  optics  provides  a  magnification  equal  to  the 
ratio  b2Vi/B1V2  where  Bi  an<*  vi  are  the  si<Jnal  bandwidths  and  acoustic 
velocity  of  the  first  cell,  etc.  This  alternate  scheme  reduces  the 
processing  rate  by  a  factor  of  about  2,  but  has  the  advantage  of 
greatly  reducing  the  bandwidth  requirement  for  the  memory  feeding  the 
second  cell.  For  a  system  composed  of  a  PbMo04  and  a  (Te02)s  Bragg 
cell,  the  memory  speed  requirement  is  reduced  from  200  Msamples/s  to 
20  Msamples/s,  as  compared  with  a  system  containing  two  PbMoOy  cells. 


3. 


IMPLEMENTATIONS. 


A  consideration  of  some  practical  limitations  and  constraints 
for  A-0  processors  and  some  experimental  measures  of  performance  were 
made  to  evaluate  to  some  degree  the  realizability  of  high  performance 
systems  based  on  the  concepts  developed  during  the  first  two  phases  of 
this  effort.  The  experimental  efforts  included  a  laboratory  bench 
version  of  a  coherent  hybrid  space/time  integrating  spectrum  analyzer 
and  a  packaged  brassboard  noncoherent  two-dimensional  ambiguity  function 
processor. 

3 . 1  _ Space/Time  Integrating  Spectrum  Analyzer. 

An  experimental  configuration  of  the  one-dimensional  coherent 
hybrid  spectrum  analyzer  of  Figure  3-1  was  implemented.  The  intent  was 
to  present  an  uncomplicated  and  straightforward  illustration  of  the 
concept  discussed  in  detail  in  the  Annual  Technical  Report  of  1980. 

The  parameters  of  the  experiment  were  chosen  to  provide  such  an 
illustration  with  devices  and  components  immediately  available  to  us, 
and  no  attempt  was  made  to  maximize  those  parameters ,  such  as  bandwidth 
and  time-bandwidth  product,  that  would  otherwise  be  important  in  the 
design  of  a  high  performance  system. 

The  optoelectronic  devices  used  in  the  experiment  were  a 
5  mW  fundamental  transverse  mode  HeNe  laser,  two  slow  shear  mode  Te02 
Bragg  cells  and  a  Reticon  CCPD1024  charge  coupled  photodiode  array. 

The  portions  of  the  bandwidth  and  time  delay  of  the  Bragg  cells  that 
were  utilized  in  the  experiment  were,  respectively,  5.3  MHz  and  6  us. 

3.1.1  _ Chirp  Spatially  Distributed  Local  Oscillators. 

One  of  the  Bragg  cells  in  this  interferometric  system  is 
driven  by  a  reference  signal  that  contains  all  the  frequencies  in  the 
analysis  band,  in  order  to  generate  a  coherent  optical  reference  wave 


3-1 


Continued 


3.1.1 


to  interfere  with  the  signal  optical  wave  at  the  detector.  It  is 
desired  that  this  reference  have  a  uniform  power  density  across  the 
band.  The  reference  waveform  used  in  these  experiments  was  a  periodic 
chirp,  generated  digitally,  with  a  variable  bandwidth  (up  to  about 
20  MHz)  and  a  variable  chirp  time.  The  chirp  time  was  set  a  6  ys. 

The  baseband  chirp  was  mixed  single  sideband  with  a  50  MHz  local 
oscillator  as  the  reference  waveform  driving  the  Bragg  cell.  The 
optical  power  spectra  of  chirps  with  bandwidths  5.3  MHz,  10.6  MHz  and 
21.2  MHz  are  shown  in  Figure  3-2  as  detected  by  the  CCD  photodiode 
array.  These  exhibit  the  familiar  Fresnel  ringing  pattern  that  is  most 
pronounced  for  small  time  bandwidth  products  (BT) .  For  comparison. 
Figure  3-3  shows  the  calculated  moduli  of  the  Fourier  transforms  of 
chirps  with  BT  values  of  20,  21,  100  and  300.  Note  that  the  details 
of  the  ringing  structure  are  very  sensitive  to  BT,  as  evidenced  by  the 
patterns  for  BT  =  20  and  BT  =  21.  The  power  spectrum  of  a  chirp  in 
practice  may  be  degraded  by  amplitude  and  phase  errors  in  the  chirp 
temporal  waveform.  These  types  of  errors  are  usually  negligible  for 
digitally  generated  chirps.  Achievable  bandwidths,  however,  are 
limited  by  the  rise-time  limitations  of  digital  systems.  A  voltage 
controlled  oscillator  (VCO)  can  produce  very  wide  band  chirps,  but 
linearity  problems  are  severe  and  chirp-to-chirp  coherence  (a  necessity 
for  these  applications)  is  extremely  difficult  to  achieve.  Surface 
acoustic  wave  (SAW)  linear  dispersive  filters  are  approaching  adequate 
amplitude  and  phase  behavior  and  are  highly  periodic  (chirp-to-chirp 
coherent) ,  controlled  by  the  timing  of  the  RF  pulses  feeding  the 
filter.  Figure  3-4a  shows  the  spectrum  (using  an  HP  spectrum  analyzer) 
of  a  RAC  SAW  filter  of  500  MHz  bandwidth  and  0.5  ys  dispersive  delay. 
Figure  3-4b  shows  a  portion  of  the  spectrum  at  an  expanded  scale;  the 
comb  structure  illustrates  the  periodic  nature  of  the  chirp  (chirps 
were  made  contiguous  with  an  0.5  ys  period). 


3-3 


Acousto-Optic  Measui 
Generated  Chirps. 


of  Power  Spectrum  of  Digitally 


Figure  3-2. 


Calculated  Modulus  of  the  Fourier  Transform  of  a  Chirp 


a)  Horizontal  Scale:  100  MHz/div. 
Vertical  Scale  :  10  dB/div. 


Figure 


b)  Horizontal  Scale:  10  MHz/div. 
Vertical  Scale  :  10  dB/div. 

3-4.  Spectrum  of  a  SAW  Dispersive  Filter. 
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3.1.2 _ Experimental  Results. 


For  the  spectrum  analyzer  experiment  a  0-5.3  MHz  chirp  was 
generated  digitally  in  6  ys  and  mixed  single  sideband  with  a  50  MHz 
local  oscillator  as  the  reference  waveform  driving  one  Bragg  cell. 

The  other  Bragg  cell  was  driven  by  a  single  tone  within  the  analysis 
band  (-50-55.3  MHz).  Figure  3-5  shows  the  integrated  intensity  in  the 
detector  plane,  as  reflected  in  the  serially  addressed  output  of  the 
photodiode  array.  Figure  3-5a  shows  the  bias  over  the  band  due  to  the 
Fourier  transformed  chirp  reference  and  the  localized  interference 
fringe  pattern  on  an  additional  localized  bias  generated  by  the  pure 
tone  signal.  In  this  experiment  BT  =  32,  so  that  the  chirp  spectrum 
nonuniformity  is  quite  pronounced.  Nevertheless,  the  spatial  frequency 
content  is  near  spatial  baseband.  Highpass  filtering  of  the  output 
shown  in  Figure  3-5a  results  in  the  signal  shown  in  Figure  3-5b,  where 
the  reference  bias  structure  and  that  due  to  the  signal  have  been 
removed. 


The  Reticon  CCPD1024  integrates  the  incident  optical  power 
density  for  250  ys,  about  the  shortest  time  in  which  the  CCD  register 
can  transfer  the  1024  charge  elements.  Since  the  Bragg  cell  delay 
and  chirp  period  are  6  ys,  we  are  integrating  by  a  factor  of  42  longer 
than  prescribed.  The  effects  of  such  a  long  exposure  is  to  introduce 
a  frequency  comb  filter.  The  teeth  of  the  comb  response  are  separated 
by  1/T  =  167  kHz  and  the  width  of  a  tooth  is  1/T^  =  4  kHz,  which  was 
observed.  When  the  frequency  is  tuned  to  a  peak,  the  intensity 
integrated  for  T^  is  (T^/T)  times  the  intensity  integrated  for  time  T. 

The  1980  Annual  Technical  Report  presented  an  analysis  of 
the  dynamic  range  characteristics  of  direct  detection  and  interfero¬ 
metric  detection,  demonstrating  the  advantages  of  the  latter.  Figure 
3-6  shows  a  comparison  of  the  measured  detector  outputs  for  direct 
detection  and  for  interferometric  detection.  For  Figure  3-6a  the 
reference  path  is  blocked  so  that  only  the  light  diffracted  from  the 
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Figure  3-5.  Signal  Plus  Reference  Spectrum  and  High  Pass 
Filtered  Version. 
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Figure  3-6.  Dynamic  Range  Comparison. 
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signal  Bragg  cell  illuminates  the  detector.  The  top  trace  corresponds 
to  a  signal  level  that  just  saturates  the  detector  and  the  lower  traces 
correspond  to  20  dB  and  30  dB  attenuation  of  the  signal.  These  traces 
were  obtained  with  two  forms  of  noise  reduction;  a  single  frame  of  the 
detector  output  in  the  absence  of  a  signal  was  subtracted  from  the 
detector  output  in  the  presence  of  the  signal,  and  the  result  was 
passed  through  a  lowpass  filter.  The  background  subtraction  served  to 
reduce  the  effects  of  diode-to-diode  nonuniformities  in  dark  current, 
and  the  lowpass  filter  constrained  the  detector  noise  bandwidth  to  the 
signal  bandwidth.  Figure  3-6b  shows  the  results  when  the  reference 
wave  was  added.  Background  subtraction  and  bandpass  filtering  were 
performed.  A  signal  is  seen  to  be  still  detectable  when  it  is  50  dB 
below  the  maximum  signal  (upper  trace)  for  interferometric  detection, 
compared  with  a  barely  detectable  signal  30  dB  down  for  direct  detection. 

The  spectral  line  shape  for  interferometric  detection  is 
shown  in  Figure  3-7.  Figure  3-7a  was  obtained  by  interfering  the  , 

signal  spectral  line  with  a  reference  wave  that  was  generated  by  driving 
the  reference  Bragg  cell  with  the  same  signal  tone  and  spreading  the 
diffracted  light  uniformly  over  detector.  The  reference  wave  has, 
therefore,  a  Doppler  frequency  shift  that  is  independent  of  position. 

The  interference  pattern  has  approximately  a  sine  envelope,  following 
the  shape  of  the  spectral  amplitude  of  the  signal,  a  result  of  the 
spatially  integrated  Fourier  transform.  This  would  also  result  from  a 
snapshot  exposure  when  the  reference  wave  is  not  the  same  frequency  as 
the  signal.  A  purely  time  integrating  spectrum  analyzer  exhibits  the 
same  spectral  line  shape,  a  result  of  integration  over  the  exposure 
time  window.  Figure  3-7a  illustrates  the  resulting  line  shape  with  the 
hybrid  space/time  integrating  technique.  The  rapidly  decaying  sidelobes 
are  evident.  This  data  demonstrates  the  characteristics  derived 
theoretically  that  the  sidelobes  for  a  coherent  hybrid  spectrum  analyzer 
are  lower  by  a  factor  of  two  in  dB  compared  with  a  purely  space- 
integrating  system,  whether  the  latter  is  interferometric  or  square-law. 
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Figure  3-7.  Spectral  Side  Sidelobe  Comparison. 
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3.2 _ Performance  Evaluation  of  a  Time  Integrating  A-0  Processor. 

A  two  dimensional  implementation  of  a  time  integrating 
acousto-optic  signal  processor,  configured  to  compute  the  cross 
ambiguity  function,  has  been  quantitatively  evaluated  in  a  limited 
series  of  tests.  Although  the  tests  are  limited  in  scope,  the  results 
provide  significant  insights  into  the  performance  limitations  and 
capabilities  of  a  real  implementation.  The  performance  results  are 
compared  with  those  of  an  ideal  cross-correlation  processor. 

3.2.1 _ Ambiguity  Function  Processor  Description. 

The  ambiguity  function  processor  is  based  on  a  time  inte¬ 
grating  acousto-optic  signal  processing  architecture  developed  at  ESL. 
Time  integrating  acousto-optic  processors  realize  flexible,  multi¬ 
purpose  complex  signal  processing  architectures  based  on  correlation 
algorithms.  Interest  in  time  integrating  techniques  is  due  to  the 
attractive  device  technology  as  well  as  flexibility  in  time  integrating 
algorithms.  An  important  consequence  of  time  integrating  techniques 
is  the  ability  to  operate  on  signals  with  very  large  time-bandwidth 
without  having  to  store  the  entire  time  history  as  a  spatial  record. 
Therefore,  an  important  class  of  two-dimentional  and  multi-channel 
algorithms  may  be  performed  without  requiring  two-dimensional  spatial 
light  modulators.  Acousto-optic  devices  and  charge  coupled  image 
sensors  are  particularly  well  suited  for  time  integrating  processor 
implementation.  Further,  these  signal  processing  techniques  may  be 
realized  with  either  coherent  or  non-coherent  optical  systems. 

The  theory  of  operation  is  summarized  in  the  following 
section.  For  more  detailed  explanation  of  time  integrating  techniques, 
refer  to  the  earlier  reports  (Time  Integrating  Optical  Signal 
Processing ,  Annual  Technical  Report,  1  August  1978-31  July  1979, 
Contract  F49620-78-C-0102  and  Interim  Report  on  Time  Integrating  Signal 
Processing ,  September  1980,  Contract  F49620-78-C-0102) . 
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3. 2. 1.1  Theory  of  Operation. 

Two-dimensional  integral  transforms,  for  which  the  kernel 
is  decomposable  in  the  proper  way,  may  be  implemented  by  time  inte¬ 
grating  optical  processing.  An  important  example  is  ambiguity 
function  processing.  Consider  the  optical  realization  shown  conceptually 
in  Figure  3-8.  The  optical  train  has  a  modulated  illumination  source, 
two  acoustic  delay  line  light  modulators,  and  a  matrix  array  of 
detectors.  This  configuration  may  be  employed  for  several  functions, 
determined  by  input  signal  and  reference  waveforms. 

The  light  diffracted  by  the  first  acousto-optic  modulator 
A01,  is  diffracted  by  the  second,  A02,  in  an  orthogonal  direction. 

Both  acoustic  signals  are  imaged  on  the  detector  plane.  The  desired 
image  plane  intensity  distribution  has  a  term  proportional  to  the 
product  of  acoustic  signals,  s.^  (t-T^)  s2  (t-x2)  ,  where  =  y/v  and 
x 2  **  x/v  are  time  variables.  The  resultant  detected  output  voltage  is 
proportional  to  the  integrated  charge.  This  structure  implements 
integrals  of  the  form: 

R(t1,t2)  =  /  sQ (t)s1(t-T1)s2 (t-x2)dt 

T 

In  this  configuration,  all  signals  are  real.  In  order  to 
perform  complex  computation,  the  real  and  imaginary  components  are  used 
to  modulate  a  carrier  in  quadrature. 


The  cross  ambiguity  function  may  be  implemented  using  the 
two-dimensional  optical  system  shown  in  Figure  3-8.  The  cross  ambiguity 
function  between  complex  signals  x(t)  and  y(t)  may  be  defined  as 


Axy(T,f)  = 


L 


x  {t)y*  (t-t  )e  j^Tfftdt 
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where  the  variables  T  and  f  may  be  interpreted  as  delay  and  Doppler, 
respectively.  The  spectral  analysis  operation  is  performed  using  the 
chirp  algorithm.  Define  the  complex  input  signals  sQ (t)  ,  s^t),  and 
s2(t)  by 

2 

-lonrt 


s0(t) 

x  (t)e 

s-^t)  = 

y*  (t) 

s2(t)  = 

eiotirt 

is  given 

by 

R<T1 ' T2 ) 

-  f 

-icnrt 

x(t)e  y*(t-x1)e 


ionr  ( t — T  2 ) 


dt 


lain . 


=  e 


» 

J  x (t)y*  (t-x1) e 


-i2iratT. 


'dt 


ICt ITT  , 


=  Axy(T1,ax2)e 


The  quadratic  phase  term  may  be  ignored  if  magnitude  square  calculation 
is  desired.  The  variables  and  correspond  to  delay  and  Doppler, 
respectively.  The  Doppler  resolution  is  commensurate  with  the  inte¬ 
gration  time;  the  analysis  bandwidth  is  determined  by  the  chirp 
bandwidth  and  acoustic  delay.  A  coarse  resolution  may  be  maintained 
during  a  signal  acquisition  period,  and  a  higher  resolution  zoom  can  be 
achieved  by  longer  integration.  Multiple  correlation  peaks  or  targets 
can  be  processed  simultaneously  since  a  linear  system  implementation  is 
used. 
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3. 2. 1.2  Acousto-Optic  Subsystem. 

The  time  integrating  signal  processor  has  been  implemented 
with  a  non-coherent  optical  system  using  a  diode  source  as  shown  in 
Figure  3-8.  Devices  include  a  Hitachi  HLP-40  light  emitting  diode, 
Fairchild  CCD221  charge  coupled  image  sensor,  and  Isomet  acousto-optic 
devices.  The  diode  has  a  30  MHz  3-dB  bandwidth  and  is  biased  at  an 
average  optical  power  of  approximately  20  mW.  The  image  sensor  has 
488  x  380  elements  and  is  operated  with  an  integration  period  of 
33.33  msec  (30  frames  per  second).  The  acousto-optic  device  have  a 
30  MHz  1-db  bandwidth  and  both  use  a  delay  range  of  approximately  20 
ysec. 


In  order  to  realize  complex  computation  with  real  input 
modulation,  the  complex  signals  modulate  a  local  oscillator  in  quadrature. 
An  interferometric  technique  with  electronic  reference  is  used.  The 
diode  is  internally  modulated  with  a  current  driver  (at  a  carrier) ,  and 
the  acoustic  devices  are  modulated  at  a  carrier  with  a  summed  reference 
oscillator.  The  analog  signal  processor  requires  three  inputs:  input 
signal  1  multiplied  by  a  chirp  waveform,  input  signal  2,  and  the 
reference  chirp  waveform.  The  required  frequency  translations  and 
chirp  multiplication  are  performed  in  an  analog  frequency  translation 
unit.  The  three  inputs  have  the  form 


Re 


Re 


Re 


i  (27Tf  t— CtTTt2) 

x  (t)  e 

i2irf .  t 

y  (t)e 

-i  (2irf  2t-onrt2)l 


where  x(t)  and  y(t)  are  the  complex  input  time  series,  fo,  f^,  and  f2 
denote  carrier  frequencies,  and  a  is  the  chirp  rate. 
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3. 2. 1.3  Bias  Elimination. 

An  optical  bias  results  due  to  the  use  of  intensity 
modulation.  Filtering  or  subtraction  may  be  used  to  eliminate  the  bias. 
The  use  of  intensity  modulation  offers  a  number  of  advantages  in 
comparison  with  coherent  systems  including  compact  non-coherent  optical 
system  realization  and  increased  immunity  to  artifact  noise.  Through 
interferometric  detection  and  complex  demodulation,  the  complex  ambiguity 
function  may  be  computed.  In  the  current  unit  a  bias  subtraction  scheme 
is  employed,  though  other  techniques  may  be  implemented.  The  bias 
subtraction  method  utilizes  two  consecutive  ambiguity  functions  where 
the  second  calculation  has  effectively  the  opposite  polarity.  That  is 
to  say,  the  first  frame  consists  of  the  bias  plus  signal  and  the  second 
consists  of  bias  minus  signal.  The  two  frames  are  subtracted  yielding 
twice  the  signal.  The  throughput  rate  is  decreased  by  a  factor  of  2. 
However,  if  filtering  were  used,  the  bias  could  be  eliminated  with  a 
single  frame.  The  bias  elimination  is  performed  digitally  in  a  high¬ 
speed  memory  (the  DeAnza  System  IP5000) . 

The  technique  that  has  been  implemented  to  provide  two 
sequential  frames  —  one  representing  bias  plus  the  real  part  of  the 
complex  ambiguity  function  and  one  representing  bias  minus  the  real 
part  of  the  complex  ambiguity  function  —  is  suboptimum.  The  technique 
provides  an  average  signal  phase  difference  of  180°  between  sequential 
frames  (thus  allowing  bias  removal  by  subtraction)  by  allowing  a 
differential  carrier  phase  drift  rate  that  is  matched  to  the  processor 
framing  rate.  Although  this  implementation  provides  the  desired  bias 
elimination,  it  uses  phase  averaging  over  180°  for  each  frame  to  achieve 
the  phase  differences  of  180°  over  sequential  frames.  Analysis  has 
shown  that  there  is  a  minimum  loss  of  3.25  dB  processing  gain  as  a 
consequence  of  this  technique.  This  loss  may  be  easily  recovered  by 
using  a  relatively  simple  phase  switching  approach.  As  a  consequence, 
effective  processing  gain  as  measured  in  the  quantitative  performance 
tests  were  at  least  3.25  dB  lower  than  they  need  to  be. 
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Quantitative  Processor  Evaluation. 


Quantitative  evaluation  of  the  A-0  processor  operating  as  a 
gated  cross-correlator  has  been  performed  by  considering  processing 
results  from  two  distinctly  different  points  of  view.  Both  of  the 
approaches  consider  the  dual  hypothesis  detection  process  associated 
with  the  detection  of  a  peak  in  the  ambiguity  plane.  The  first 
approach  is  based  on  the  assumption  that  the  underlying  statistics 
of  the  detection  random  variables  are  unknown  and  it  is  necessary  to 
develop  the  comparison  with  an  ideal  processor  through  measured  prob¬ 
ability  of  detection,  P^,  and  probability  of  false  alarm,  Pfa> 
characteristics.  In  the  second  approach,  the  underlying  statistics  of 
the  detection  process  are  hypothesized  (and  tested  to  a  limited  extent) 
and  measured  processor  output  SNR  values  are  compared  with  those  of  an 
ideal  correlator. 

3. 2. 2.1  Test  Conditions. 

The  quantitative  tests  on  the  A-0  processor  all  used 
simulated  signals  with  additive  Gaussian  noise.  The  signals  of  interest 
have  been  simulated  as  a  repetitive  pulse  train  with  a  multipath  delay 
of  nominally  25  ms  and  with  PRI  variation  between  19  ys  and  312  ys. 

The  pulse  width  has  been  fixed  at  1  ys.  The  input  signal-to-noise 
ratios  (both  direct  signal  and  multipath)  are  specified  as  the  ratio  of 
peak  signal  power  (power  measured  during  the  1  ys  pulse)  to  the  equiva¬ 
lent  noise  power  measured  in  a  1-MHz  bandwidth.  The  pre-D  bandwidth 
for  the  composite  signal  plus  noise  was  varied  from  300  kHz  to  1  MHz. 

Gated  correlation  processing  was  used  on  all  of  the  tests. 

The  gate  was  derived  from  the  noisy  signal  in  much  the  same  way  in 
which  the  process  might  be  automated  in  a  real  multipath  processing 
configuration.  Specifically  the  noisy  pre-D  signal  is  AM  detected  and 
passed  to  a  recursive  synchronous  video  integrator  (SVI)  that  has  a 
length  corresponding  to  the  PRI  of  the  desired  signal  (assumed  known 
a  'priori )  . 
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3.2.2. 1  —  Continued 

The  output  of  the  SVI  is  reconstructed  through  a  signal 
synchronizer  to  provide  a  repetitive  gate  signal  that  is  centered  on 
the  direct  signal  and  has  a  width  on  the  order  of  2  ys.  The  width  of 
the  gate  may  be  varied  to  approximately  match  the  signal  that  has 
passed  through  a  pre-D  filter.  The  gate  signal  as  well  as  the  pre-D 
input  signal  that  is  used  to  modulate  the  LED  in  the  ambiguity  function 
processor  are  digitally  delayed  by  a  pre-computational  delay  unit  that, 
by  adjustment,  provides  the  means  to  select  the  delay  range  to  be 
processed  in  the  ambiguity  function  processor.  That  is,  the  delay 
window  of  approximately  20  ys  can  be  centered  on  any  particular  delay 
value  by  adjustment  of  the  pre-computational  delay. 

A  block  diagram  of  the  experimental  configuration  is  shown 
in  Figure  3-9.  The  ambiguity  function  processor  in  this  case  includes 
a  digital  chirp  generator,  a  frequency  translation  unit,  the  A-0  triple 
product  processor,  a  camera  control  unit,  and  a  special  purpose  digital 
image  processor.  The  Doppler  window  is  centered  on  the  desired  range 
(positive  or  negative)  by  simply  tuning  the  frequency  synthesizer.  The 
range  of  Doppler  that  is  processed  is  selected  by  choosing  the  chirp 
rate.  For  the  tests  that  are  reported  here,  the  Doppler  range  is 
approximately  1.1  kHz. 

3. 2. 2. 2  Data  Collection. 

For  each  set  of  conditions  a  block  of  data  is  collected  so 
that  subsequent  analysis  can  provide  as  complete  a  picture  as  possible 
of  the  processor  performance.  The  processor  is  set  up  with  controlled 
and  known  conditions.  These  include  signal  properties  such  as  RFI , 
pulse  width,  direct  path  and  multipath  SNR’s.  Processor  parameters 
such  as  pre-D  filtering  bandwidth,  gate  width,  Doppler  and  delay  offset, 
and  Doppler  range  are  also  established.  For  each  of  these  fixed  con¬ 
ditions,  a  large  sequence  of  processing  runs  (generally  on  the  order  of 
1000,  however  some  at  10,000)  are  performed  to  measure  the  statistical 
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Figure  3-9. 


Experimental  Configuration  for  Repetitive  Pulse  Signal 
Multipath  Processing. 
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3. 2. 2. 2  —  Continued 


properties  of  the  ambiguity  function  envelope  in  the  noise  region  and 

in  the  signal  peak  region.  The  properties  include  the  mean,  the 

standard  deviation  and  the  cumulative  distribution  function.  The  number 

of  points  used  for  the  noise  only,  and  signal  plus  noise  CDFs  were  10,000 

and  1,000  respectively.  At  a  selected  threshold  value,  the  P_  is  given 

by  one  minus  the  noise  only  CDF  value  and  is  given  by  one  minus  the 

signal  plus  noise  value.  As  an  example  from  the  table,  a  threshold 

_  2 

number  of  108  corresponds  to  a  P^a  of  10  and  a  Pd  of  0.39.  The 
additional  statistics,  mean  and  standard  deviation,  are  useful  in  the 
performance  analysis  in  a  way  that  will  be  described. 

3.2.2. 3  Processor  Analysis  Based  on  P,  and  P . . 


As  has  been  described  in  the  preceeding  section,  the  CDFs  of 
the  linearly  detected  envelope  at  the  peak  of  the  cross  ambiguity 
function  and  in  a  "noise  only"  region  of  the  ambiguity  function  plane 
have  been  measured  for  various  input  signal  conditions.  By  varying 
the  threshold  on  the  pair  of  CDFs  a  curve  of  Pd  versus  Pfa  may  be  con¬ 
structed.  This  type  of  curve  is  normally  called  a  ROC  curve  or  Receiver 
Operating  Characteristic. 

Theoretical  ROC  curves  for  a  detection  process  based  on  a 
decision  at  the  output  of  a  linear  envelope  detector  where  the  input 
is  either  Gaussian  noise  only  or  a  sinusoid  in  Gaussian  noise  are 
available  in  the  literature.  Curves  of  this  type  will  enable  us  to 
relate  the  measured  P^/P^  characteristics  to  an  equivalent  input  SNR 
for  the  detection  of  a  single  sinewave  pulse  in  Gaussian  noise. 

Receiver  operating  characteristics  for  which  values  of 
constant  SNR  are  straight  lines  have  been  developed  by  Robertson  (G.H. 
Robertson,  "Operating  Characteristics  for  a  Linear  Detector  of  CW 
Signals  in  Narrow-Band  Gaussian  Noise",  BSTJ,  April  1967,  pp  755-774). 
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An  example  of  the  set  of  curves  is  given  in  Figure  3-10  (from  A.D. 

Whalen,  Detection  of  Signals  in  Noise,  Academic  Press,  Mew  York,  1971, 
p.  248).  Comparable  data  are  also  available  in  D.K.  Barton,  Radars , 
Volume  2,  The  Radar  Equation,  Artech  House,  Inc.,  1974,  p.  115. 

The  theoretical  results  represented  by  these  curves  are 
based  on  the  detection  of  the  envelope  of  a  single  pulse  of  sinusoid  in 
Gaussian  noise.  Thus,  as  we  superimpose  the  measured  ROC  data  points, 
departure  of  the  set  of  points  from  a  line  parallel  to  those  for  con¬ 
stant  input  SNR  could  be  indicative  of  departure  from  the  Gaussian 
assumption.  Assuming  there  is  nominally  a  straight  line,  the 
"equivalent  single  pulse  input  SNR",  Yq  maY  be  estimated  by  curve 
matching  and  interpolation.  Measured  data  values  for  a  variety  of  input 
conditions  are  shown  superimposed  on  the  theoretical  ROC  curves  in 
Figure  3-11. 

These  data  may  be  compared  with  an  "equivalent  single  phase 
input  SNR",  rQ,  that  has  been  derived  from  an  analysis  of  an  ideal 
gated  correlator.  For  idealized  gated  correlation  processing, 

(SNRX) (SNR2) 

F0  =  n  (SNR^  +  1) 

where  SNR^  is  the  direct  path  input  SNR,  SNR2  is  the  multipath  input 
SNR,  gating  is  performed  on  the  direct  path  pulses  and  n  of  these  pulses 
are  used  in  the  processing. 

The  results  of  the  data  analysis  are  summarized  in  Figure  3-12. 
Based  on  this  comparison,  the  apparent  loss  in  processing  gain  relative 
to  an  ideal  gated  cross-correlator  ranges  from  10.6  dB  to  16.8  dB.  For 
these  data,  the  largest  apparent  loss  occurs  for  conditions  under  which 
the  expected  processing  gain  is  largest.  The  possible  increase  in 


Gated  Correlator 


Measured  SNR,  yQ  (dB) 


A  Comparison  of  Processor  Equivalent  SNR  for  an  Ideal 
Gated  Cross-Correlator  (Parker's  Model)  and  Corresponding 
Values  Derived  from  Measured  P./Pf  Characteristics  of 
the  A-0  Processor. 


Figure  3-12. 
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3. 2. 2. 3 


apparent  loss  for  high  processing  gain  (a  larger  number  of  pulses  in 
the  processing  interval)  will  be  examined  further  in  a  subsequent 
section. 


Another  possible  insight  into  processor  performance  may  be 
available  by  considering  the  apparent  loss  as  a  function  of  the  multi- 
path  SNR.  These  data  are  shown  in  Figure  3-13.  A  straight  line  with 
a  slope  of  1/2  is  superimposed  on  the  data  to  illustrate  a  trend. 


3. 2. 2. 4  Detector  Input  Output  Relationships. 


The  pd/pfa  measurements  have  been  used  to  provide  an  estimate 
of  the  effective  SNR,  -y^,  at  the  input  to  a  linear  envelope  detector. 

As  part  of  our  data  collection,  the  mean  and  variance  of  the  output  of 
the  linear  envelope  detector  have  also  been  measured  both  in  a  noise 
region  and  at  the  peak  of  the  ambiguity  function.  Based  on  the  central 
limit  theorem  we  might  expect  the  random  part  of  our  ambiguity  function 
estimate  (before  taking  the  envelope  and  in  any  region  of  the  plane)  is 
Gaussian.  The  fact  that  the  measured  ROC  data  points  generally  follow 
straight  lines  (Figure  3-11)  supports  a  Gaussian  assumption.  Additional 
support  is  provided  by  the  measured  statistics.  Specifically,  the 
envelope  of  a  zero-mean  Gaussian  random  variable,  as  is  expected  in  the 
"noise-only"  region  of  the  ambiguity  plane,  is  Rayleigh  and  the  ratio 
of  the  mean  to  the  standard  deviation  is 
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For  a  set  of  21  measurements  (each  of  which  represents  the  computation 
of  a  thousand  or  more  ambiguity  surfaces) ,  the  mean  of  this  ratio  is 
1.904  and  the  standard  deviation  is  0.047. 
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Figure  3-13.  Apparent  Loss  as  a  Function  of  Multipath  Input  SNR  (dB) 
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For  a  linear  envelope  detector  with  a  sinewave  plus  zero 
mean  Gaussian  noise  on  the  input,  the  first  and  second  moments  of  the 
output,  z,  are  given  by  Whalen  (Anthony  D.  Whalen,  Detection  of 
Signals  in  Noise,  Academic  Press,  New  York,  1971,  p  105)  as 

2 

£{z}  =  (2o2 )  (j)  (— 1;  -^jr) 

2 

£{z2}  =  (2a2)T  (2)1F1(-1;  1,  -5j-) 

2 

where  o  is  the  variance  of  the  Gaussian  noise,  a  is  the  ratio  of  the 
sinewave  amplitude  to  the  rms  level  of  the  noise,  TO  is  a  gamma 
function  and  b;  x)  is  the  confluent  hypergeometric  function. 

If  one  computes  the  detector  output  SNR  as 


where  u  is  the  mean  detector  output  and  cj_  is  the  standard  deviation 
s  s 

of  the  detector  output  (both  taken  at  the  peak  of  the  ambiguity 
function),  then  at  high  SNR  values  (say  greater  than  10  dB,  the  output 
SNR  equals  the  input  SNR.  The  complete  relationship  is  shown  in 
Figure  3-14.  Superimposed  on  the  theoretical  curve  are  measured  data 
developed  from  the  output  envelope  statistics  and  Yq  (developed  from 
the  Pd/Pfa  measurements) .  These  data  provide  additional  support  for 
the  contention  that  to  establish  performance  we  can  simply  measure 
statistics  at  the  peak  of  the  ambiguity  function  and  from  these  infer 
expected  Pcj/Pfa  characteristics. 
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3. 2. 2. 5  Processor  Evaluation  Based  on  Peak  Envelope  Statistics. 


For  each  of  the  various  measurement  conditions ,  the  measured 
mean  and  standard  deviation  values  for  the  envelope  of  the  peak  of  the 
ambiguity  function  can  be  used  to  compute  an  effective  output  SNR  that 
in  turn  can  be  transformed  through  the  theoretical  curve  of  Figure  3-14 
to  an  equivalent  input  SNR  that  will  be  called  y^.  These  data  will  be 
compared  with  theoretical  results  based  on  an  analysis  of  an  ideal 
cross-correlation  processor.  The  model  in  this  case  was  reported  in 
the  Interim  Technical  Report,  1979.  The  ideal  processor  output  SNR  is 
defined  as  the  ratio  of  the  squared  mean  of  the  correlation  estimate  to 
its  variance.  The  output  SNR,  r ^ ,  is  given  by 

2B  T  SNR.  SNR., 
p  — .  n  X  ^ 

1  ~  1  +  SNRX  +  SNR2  +  XSNR1SNR2 

where  2BnT  is  the  time-bandwidth  product  and  is  in  our  case  equal  to 
n,  the  number  of  pulses;  X  is  a  parameter,  typically  with  a  value  near 
unity.  The  X  term  takes  into  account  the  estimation  error  due  to 
finite  time  integration. 
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The  measured  results  have  been  compared  with  those  of  an 
ideal  correlator  for  values  of  X  equal  to  one  or  zero  and  are  presented 
in  Figure  3-15.  The  data  seem  to  be  concentrated  in  a  band  with  an 
apparent  loss  ranging  from  10  to  13  dB  and  a  second  band  with  an 
apparent  loss  on  the  order  of  16  dB.  If  the  data  are  examined  in  terms 
of  the  apparent  loss  as  a  function  of  n,  the  results  are  summarized  in 
Figure  3-16.  Here  it  seems  clear  that  the  larger  apparent  losses  are 
occurring  when  n  is  large.  Large  n  corresponds  to  the  higher  duty 
cycle  signals,  in  this  case  a  PRI  of  about  19  ysec.  Note  that  the  number 
of  pulses  also  corresponds  to  the  processing  time-bandwidth  product  or 
equivalently  the  processing  gain.  Performance  analysis  of  the  optical 
correlator,  as  reported  in  the  1979  Interim  Report  has  suggested  a 
processing  gain  of  the  form 
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n 


1  +  k 


n 


(DR) 


2 


where  k  is  a  constant  related  to  the  signal  modulation  depths  of  the 
three  optical  processor  inputs  and  DR  is  the  dynamic  range  of  the  photo¬ 
detector  array.  For  very  large  n  the  processing  gain  approaches  a 

2 

constant  given  by  (DR)  /k.  At  this  time,  we  have  insufficient  measure¬ 
ments  to  establish  the  presence  of  the  effect,  however,  it  is  a 
candidate  for  controlled  experiments  in  the  future. 

The  apparent  loss  is  plotted  as  a  function  of  the  multipath 
input  SNR  in  Figure  3-17.  These  data  may  be  further  compared  with  the 
data  presented  in  Figure  3-13.  With  the  additional  data  available  in 
Figure  3-17  the  possible  trend  identified  in  Figure  3-13  is  not  apparent. 

3, 2. 2. 6  Summary  of  Quantitative  Test  Results. 

The  time  integrating  triple  product  processor  has  been 
configured  as  an  ambiguity  function  processor  and  operated  in  a  gated 
processing  mode  to  evaluate  its  performance  for  the  detection  of 
multipath.  The  test  signal  is  a  repetitive  pulse  train  with  a  1  ysec 
pulse  width  and  a  PRI  that  was  varied  between  19  ysec  and  312  ysec. 

The  multipath  delay  was  nominally  25  ysec  and  the  pre-D  bandwidth  for 
the  composite  signal  plus  noise  was  varied  from  300  kHz  to  1  MHz. 

Gated  correlation  processing  allows  one  to  increase  the 
effective  average  SNR  for  pulse  signals.  The  gate  in  this  case  was 
derived  from  the  noisy  pre-D  signal  and  had  a  nominal  width  of  2  ysec. 

A  key  attribute  of  the  processor  is  its  highly  parallel 
processing  capability.  For  the  tests  reported  here  the  processor  was 
configured  to  compute  the  real  part  of  the  complex  cross  ambiguity 
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Apparent  Loss  as  a  Function  of  Multipath  Input  SNR  (db) 
X  =  1  Model. 
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function  in  real  time  and  for  a  region  20  ysec  by  1.1  kHz.  The 
measurement  interval  was  (1/15)  sec,  however,  the  effective  integration 
time  depended  upon  the  gate  width  and  the  signal  PRI .  The- region  of 
the  ambiguity  function  that  is  processed  may  be  centered  wherever 
desired  or  sequentially  stepped  in  a  search  scenario  by  the  simple 
selection  of  a  pre-computational  delay  parameter  value  and  a  synthesizer 
frequency  parameter  value. 

The  processor  output  values  are  in  digital  form  and  at  a 
cursor  selectable  position  a  series  of  digital  samples  are  used  to 
compute  the  envelope  of  the  complex  ambiguity  function  via  the  Hilbert 
transform. 


The  performance  evaluation  was  placed  in  the  context  of  a 
detection  problem,  that  is,  is  signal  plus  noise  present  or  is  noise 
only  present.  The  signal  here  is  the  real  part  of  the  complex 
ambiguity  function  at  a  carrier  and  at  the  multipath  delay  and  Doppler 
offset.  The  processor  P^/P^  characteristics  were  estimated  by 
performing  a  large  series  of  tests  at  various  input  SNR's.  Additional 
statistical  characteristics  such  as  the  mean  and  variance  of  the 
envelope  output  were  measured  in  both  signal  and  noise  only  regiois. 

The  measured  data  were  processed  to  arrive  at  an  effective 
SNR  of  a  sinewave  (phase  unknown)  in  Gaussian  noise  that  would  be 
applied  to  a  linear  envelope  detector  for  the  detection  process.  This 
effective  SNR  is  a  consequence  of  the  processor  and  can  be  compared 
with  an  equivalent  SNR  for  an  ideal  correlator  operating  with  the  same 
input  signals.  The  differences  between  the  effective  SNR's  derived 
from  measurements  and  comparable  values  for  the  ideal  correlator  repre¬ 
sent  an  apparent  loss  in  processing  gain  associated  with  this  particular 
real  processor  implementation. 
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The  apparent  losses  that  were  measured  concentrated  in  a 
band  ranging  from  10  to  13  dB  and  in  a  second  band  with  an  apparent 
loss  on  the  order  of  16  dB.  The  higher  apparent  loss  was  observed  at 
high  processing  time-bandwidth  product  and  may  be  a  manifestation  of 
an  effect  anticipated  by  Kellraan  and  associated  with  a  signal  modulation 
depths  in  the  optical  processor  and  dynamic  range  limtiations  in  the 
photo-detector  array.  It  has  been  determined  that  a  minimum  of  3.25  dB 
of  this  apparent  loss  is  due  to  the  phase  averaging  technique  used  for 
bias  elimination.  Alternative  bias  elimination  techniques  are  available 
that  do  not  suffer  this  loss  and  that  are  relatively  simple  to  implement. 
With  the  assumption  that  this  is  correct  and  if  we  consider  the  low-duty 
cycle  signals  (lower  effective  processing  gain  available  in  the  measure¬ 
ment  interval) ,  then  it  is  concluded  that  the  apparent  loss  of  the 
optical  processor  relative  to  ideal  is  on  the  order  of  7  to  10  dB.  The 
level  of  processing  gain  that  is  achievable  is  quite  satisfactory  for 
many  applications.  If  the  apparent  loss  could  be  reduced  by  another 
3  to  6  dB  many  more  applications  would  become  viable.  It  is  recommended 
that  additional  controlled  tests  be  performed  to  establish  the  sources 
of  the  loss  and  evolve  techniques  to  improve  performance. 
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DETECTOR  INTEGRATION  ACOUSTO-OPTIC 
SIGNAL  PROCESSING 

P.  KELLMAN 
ESL  Incorporated 
Sunnyvale,  California  94086 


1.  INTRODUCTION 

Detector  integration  optical  signal  processing 
techniques  presented  in  this  paper  are  a  generali¬ 
zation  of  the  one  dimensional  time  integrating  cor¬ 
relator  demonstrated  by  Sprague1  and  transversal 
filter  architectures  described  by  Whitehouse, 
et  al2.  The  class  of  achievable  signal  processing 
functions  is  described,  and  examples  of  one  and  two 
dimensional  realizations  of  spectral  analysis  and 
ambiguity  function  processing  are  given. 

Detector  integration  optical  processing  methods 
implement  time  integrating,  rather  than  spatial 
integrating,  algorithms,  which  result  in  a  processor 
with  increased  flexibility  and  multiple  purpose.  An 
important  consequence  of  time  integrating  methods  is 
the  ability  to  operate  on  signals  or  imagery  with 
very  large  time-bandwidth  products  (M.06 )  without 
having  to  store  the  entire  time  history  as  a  spatial 
record.  Therefore,  an  important  class  of  two- 
dimensional  and  multi-channel  processing  algorithms 
may  be  performed  without  requiring  two-dimensional 
spatial  light  modulator  devices.  Rather,  the  opti¬ 
cal  implementation  presented  utilizes  acousto-optic 
devices  and  charge  coupled  image  sensors,  which  are 
highly  compatible  input/output  devices  for  signal 
processing  applications.  Further,  these  techniques 
nay  be  realized  with  either  coherent  or  non¬ 
coherent  optical  systems. 

The  one  dimensional  time  integrating  correla¬ 
tor  is  reviewed  in  Section  2  and  a  spectral  analy¬ 
sis  implementation  is  described  in  Section  3.  Two 
dimensional  detector  integration  processors  are 
introduced  in  Section  4  and  examples  of  ambiguity 
function  processing  and  two  dimensional  spectral 
analysis  are  given  in  Sections  5  and  6. 

2.  TIME  INTEGRATING  CORRELATOR 

The  one  dimensional  time  integrating  correla¬ 
tor1  is  reviewed.  Consider  the  optical  realization 
shown  conceptually  in  Figure  1.  The  light  source 
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is  temporally  modulated  by  a  signal,  s^lt),  to 
produce  am  output  intensity  I^(t),  which  illumi¬ 
nates  an  acoustic  delay  line  light  modulator.  A 
second  signal,  s2<t),  is  introduced  into  the 

acoustic  cell,  which  spatially  modulates  the  source 
intensity  by  an  amount  denoted  by  I2(t-x/v),  where 

x  is  the  spatial  dimension  and  v  is  the  acoustic 
velocity.  The  acoustic  signal  is  Schlieren  imaged 
onto  a  linear  photodiode  array;  the  intensity  dis¬ 
tribution  in  the  image  plane  is  given  by  the  pro¬ 
duct  (t) I2 (t-x/v) . 

The  charge  integration,  directly  proportional 
to  exposure,  is  performed  by  detectors  at  discrete 
positions,  xi*  The  resultant  output  voltage  at  the 

ith  detector  element  is  therefore; 

R12(V  “  A  Vm2<t''xi/V)dt  U) 

where  -  Xj/v,  and  T,  the  integration  time  of  the 

detector,  is  set  by  the  timing  of  the  charge  trans¬ 
fer  readout  register.  This  is  the  desired  correla¬ 
tion,  where  and  are  directly  related  to  the 

input  signals  s^  and  s2> 

The  range  of  relative  delay  between  input 
signals  is  limited  by  the  acoustic  delay  length, 
and  the  integration  time  is  set  by  the  detector 
readout  period.  In  the  spatial  integration 
approach1,  the  integration  period  is  determined  by 
the  acoustic  delay;  in  configurations  employing 
fixed  reference  masks,  the  range  of  relative  delay 
is  unlimited,  however,  in  correlators  with  active 
reference,  the  range  is  limited  to  the  acoustic 
delay.  Long  integration  and  flexibility  of  vari¬ 
able  integration  are  achieved  with  the  time  inte¬ 
grating  approach. 

In  order  to  operate  in  a  linear  region  of 
intensity  modulation,  it  is  necessary  to  bias  the 
input  signals.  The  resulting  intensity  modulation 

is; 


I2  =>  sin2  [B  +  Sjl  =<  B2  +  s2  (2) 

where  the  acousto-optic  modulation  is  considered 
in  the  Bragg  limit.  In  this  case  the  spatial  modu¬ 
lation  is  approximately  linear  for  weak  modulation 
biased  at  50*  diffraction  efficiency.  The  output 
voltage  may  then  be  written  as, 


Figure  1.  Time  Integrating  Correlator 


R,  <  T  .  )  =  fl,  <t)I, (t-x./v)dt 

12  x  1  2  1 

=  ^1(t)s2(t-xi/v)dt  +  ./BjBjdt 

+  B^  t-x^/v) dt  +  B2  /s1(t)dt  .  (3) 

The  first  term  of  Equation  (3)  is  the  desired  cor¬ 
relation  between  and  s^.  The  second  term  is 

fixed  bias,  and  the  last  terms  are  signal  depend¬ 
ent,  though  effectively  zero  mean  (signal  com¬ 
ponents  at  frequencies  higher  than  1/T  will 
integrate  to  zero) . 

Operation  at  low  signal  to  bias  ratio  (modu¬ 
lation  depth) ,  in  order  to  achieve  linear  intensity 
modulation,  reduces  the  correlator  dynamic  range. 
Furthermore,  wideband  acousto-optic  devices  with 
high  diffraction  efficiency  are  difficult  to 
realize.  A  technique,  introduced  in  this  paper, 
for  overcoming  this  limitation,  is  described  as 
follows.  The  acousto-optic  deflector  is  modulated 
single  sideband  with  the  carrier  offset  by  an 
amount  greater  than  the  signal  bandwidth;  in  addi¬ 
tion  the  illumination  source  must  be  modulated  on 
a  carrier  with  equal  frequency  offset.  In  this 
way,  undesirable  products  resulting  from  non¬ 
linear  intensity  modulation  are  offset  in  frequency, 
and  integrated  to  zero.  Consequently,  low  signal 
to  bias  ratio  and  high  diffraction  efficiency  are 
not  required.  The  spatial  sampling  in  the  detec¬ 
tion  plane  must  then  be  commensurate  with  the  band¬ 
width  of  signal  plus  frequency  offset.  The  spatial 
carrier  can  be  eliminated  by  means  of  a  grating  in 
the  image  plane,  thereby  reducing  the  detector 
sampling  requirement.  Other,  more  general,  func¬ 
tions  can  be  realized  with  spatial  masks,  and  by 
post  detection  electronic  weighting. 

Detector  noise  determines  the  achievable  time- 
bandwidth  product  for  a  given  sensitivity.  Detec¬ 
tor  noise  has  a  shot  noise  contribution  due  to  the 
optical  bias,  as  well  as  post  integration  readout 
noise  due  to  reset  charge  uncertainty  in  charge 
transfer  devices.  The  achievable  time -bandwidth 
and  dynamic  range  can  be  extended,  however,  by 
employing  post  detection  digital  integration. 
Baseline  variation  due  to  optical  bias  and  detec¬ 
tor  array  non- uniformities  also  reduces  the  dynamic 
range.  However,  this  baseline  may  be  effectively 
removed  by  electronic  techniques. 

3.  TIME  INTEGRATING  SPECTRAL  ANALYSIS 

An  approach  to  real  time  optical  spectral 
analysis  of  electrical  signals,  that  uses  a  time 
integrating,  rather  than  spatial  integrating,  ver¬ 
sion  of  the  chirp-z  transform,  is  described.  The 
large  time -bandwidth  correlation  is  accomplished 
by  means  of  the  correlator  described  in  the  pre¬ 
vious  section.  A  time  integrating  architecture, 
for  correlation  and  spectral  analysis,  using 
oppositely  traveling  acoustic  waves,  has  been 
described  by  Montgomery1'.  For  time  integrating 
realizations,  frequency  resolution  is  determined 
by  the  detector  integration  period,  rather  than 
the  acoustic  delay  length,  therefore,  higher 
resolution  can  be  achieved  than  by  spatial  inte¬ 
gration.  Variable  resolution  is  attainable  through 


variable  detector  integration.  Use  of  non-coherent 
intensity  modulation  leads  to  detection  of  the  mag¬ 
nitude  spectrum  rather  than  the  power  spectrum 
which  yields  a  considerable  increase  in  dynamic 
range.  This  approach  to  spectral  analysis  achieves 
a  flexibility  not  readily  achieved  by  coherent 
optical  spatial  methods.  The  chirp  algorithm  is 
reviewed,  and  the  real  implementation  is  discussed. 

The  Fourier  transform  integral  of  Equation  (4) 
may  be  rewritten  as  Equation  (5)  by  expressing  f«t 

as  1/2 ( f 2  +  t2  -  (t-f)2]. 


S(f> 

S(f) 


/»  -i2TTft 
■  Js  (t)e 


-inf' 


dt 


2  2  2 
r  -ixt  in (t-f) 
ys(t)e  e  dt 


(4) 

(5) 


The  conqplex  realization  is  shown  in  Figure  2.  The 
post  phase  weighting  may  be  ignored  if  only  the 
magnitude  or  power  spectrum  is  required.  This 
realization  converts  the  integral  transform  to 
time  invariant  filtering  of  a  preweighted  signal, 
followed  by  post  weighting. 
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Figure  2.  Complex  Realization  of  Chirp  Algorithm 

In  general,  linear  time-variant  operations  of 
the  form: 

y(T)  - /x(t)h(t,T)dt  (6) 

may  be  implemented  whenever  the  kernel  can  be 
decomposed  as: 

h(t,T)  -  a(t)b(t-T)c(T)  (7) 

A  necessary  and  sufficient  condition  for  the  above 
decomposition5  is  given  as: 

Ap.9  »  +  h  .  0  (8) 

St  3t  StST 

The  chirp  algorithm  may  be  realized  with  a 
single  real  correlation  by  translating  the  real 
chirp  to  a  carrier  frequency,  f  ,  greater  than  the 

signal  analysis  bandwidth,  B.  This  implementation 
is  shown  in  Figure  3.  The  desired  difference  fre¬ 
quency  is  translated  to  baseband,  and  the  sum  fre¬ 
quency  is  an  out  of  band  chirp  that  integrates  to 
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Figure  3.  Single  Channel  Real  Implementation 
of  Chirp  Algorithm 


zero.  The  product  of  the  difference  frequency, 
which  is  a  linear  function  of  delay,  and  the  input 
signal,  is  integrated,  producing  the  spectrum. 

Frequency  resolution  of  this  spectrum  analyzer 
is  inversely  proportional  to  the  detector  integra- 
tion  time,  T.  The  resolution  is  not  set  by  the 
acoustic  time  delay  which  typically  limits  the  use 
of  acousto-optical  processing  to  wideband  analysis. 
A  wide  range  of  integration  times  are  possible 
(0.1  -  10  msec)  using  self-scanned  arrays  of 
moderate  size  (e.g. ,  1024).  Integration  time 
periods  can  be  extended  easily  by  means  of  elec¬ 
tronic  accumulation. 

The  delay  line  length,  T,  and  the  number  of 
detectors,  N,  set  the  difference.  At,  between 
detector  elements  in  the  image  plane.  This  deter¬ 
mines  the  chirp  bandwidth,  B' .  The  required  reso¬ 
lution  or  integration  period  fixes  the  chirp  rate 
(see  Figure  4).  The  above  relations  may  be  written 
as : 


T  -  N.Al 
Af  =  1/T 
B*  Af 

-  =i  - —  -*■  B 

T  At 
B'  =  N/T 


(9) 


1/At 


The  chirp  bandwidth,  B' ,  will  be  greater  than  the 
analysis  bandwidth,  B,  if  the  integration  time,  T, 
is  greater  than  the  delay  line  length,  T. 
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in  dynamic  range  as  compared  to  coherent  optical 
spectral  analysis  in  which  incident  exposure  is 
related  to  spectral  power  density  and  the  dynamic 
range  is  given  by  10  log^M. 

4.  TWO-DIMENSIONAL  PROCESSING 

Two-dimensional  integral  transforms,  for  which 
the  kernel  is  decomposable  in  the  proper  way  may  be 
implemented  by  detector  integration  optical  pro¬ 
cessing.  Two  important  exaoples  are  ambiguity 
plane  processing  and  spectral  analysis  which  are 
described  in  the  next  sections.  The  general  class 
of  achievable  operations  is  described  in  this  sec¬ 
tion. 

Consider  the  optical  realization  shown  con¬ 
ceptually  in  Figure  5.  The  optical  train  has  a 
modulated  illumination  source,  two  acoustic  delay 
line  light  modulators,  and  a  matrix  array  of  detec¬ 
tors.  This  configuration  may  be  employed  for 
several  functions,  determined  by  input  signal  and 
reference  waveforms.  It  is  assumed  throughout  this 
discussion,  that  two-dimensional  processing  is 
applied  to  either  very  long  one-dimensional  signals 
or  to  two-dimensional  signals  (e.g.,  imagery)  that 
are  in  raster  format.  In  this  way,  the  input, 
x(t) ,  will  be  a  function  of  one  variable,  t;  the 
output,  yd^.Tj),  is  a  function  of  two  variables. 
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Figure  4.  Instantaneous  Frequency  of  Chirp 

The  samples  must  be  spaced  less  than  1/T 
(modifying  Equation  9)  in  order  to  satisfy  Nyquist 
sampling  of  the  signal  plus  the  carrier  offset. 

The  spatial  carrier  may  be  eliminated  by  means  of 
a  grating  in  the  image  plane,  or  by  electronic 
detection.  The  magnitude  spectrum  corresponds  to 
the  detected  envelope)  the  phase  modulates  the 
carrier  and  may  be  detected  as  well,  resulting  in 
a  complex  spectrum  output. 

The  dynamic  range  is  primarily  limited  by  the 
detectors,  provided  there  is  sufficient  signal 
light  power.  If  there  are  M  distinguishable  levels 
in  the  output  voltage  between  the  noise  floor  and 
the  photoelement  saturation  value,  the  dynamic 
range  in  dB  is  given  by  20  log^QM,  since  incident 

exposure  is  directly  proportional  to  spectral  mag¬ 
nitude.  This  results  in  a  considerable  increase 


Figure  5.  Generalized  2-D  Detector  Integration 
Processor 

The  light  transmitted  by  the  first  acousto¬ 
optic  modulator  A01,  is  diffracted  by  the  second, 
A02,  in  an  orthogonal  direction.  Both  acoustic 
signals  are  imaged  on  the  detector  plane.  The 
desired  image  plane  intensity  distribution  has  a 
term  proportional  to  the  product  of  acoustic  sig¬ 
nals,  b(t-T^) c (t-T^) ,  where  »  x/v  and  T2  »  y/v 

are  time  variables.  The  resultant  detected  output 
voltage  is  proportional  to  the  integrated  charge  as 
in  the  previously  described  correlator.  This 
structure  implements  integrals  of  the  form: 

y(T1,T2)  -  Xf  x(t)h(T1(T2,t)dt  (10) 

where 

h(Tl,T2)  -  a(t)b(t-T2)c(t-T1)d(T1,T2)  ,  (11) 


s2<t-T1)dt 


[a,  b,  c ,  and  d}  are  reference  signals,  and  x(t)  is 
the  input  signal.  In  this  configuration,  all 
signals  are  real.  In  order  to  perform  complex  com¬ 
putation,  similar  considerations  apply  that  were 
given  for  one-dimensional  processing.  The  output 
weighting,  d(x^,l2)#  may  be  applied  electronically 

as  indicated,  or  by  means  of  a  spatial  mask  in  the 
image  plane. 

Alternate  arrangements  include  applying  inputs 
at  A01  or  A02.  A  dual  input  processor,  that  uses 
the  same  optical  system,  implements  integrals  of 
the  form: 

y(tj.t2)  =  /x1(t)x2(t-T1)h(T1,T2,t)dt  (12) 

where  h  is  given  by  Equation  (11) ,  the  source  modu¬ 
lation  is  x.(t) ,  and  A02  is  driven  by  c(t)x2(t). 

Additional  input  sources  and  multichannel  acoustic 
devices  may  be  employed  as  well.  These  possibili¬ 
ties  are  not  detailed  here. 

The  two-dimensional  time  integrating  optical 
processor  may  be  implemented  with  several  modula- 
tion/detection  schemes.  A  description  of  a  coher¬ 
ent  optical  implementation  using  an  interferometric 
reference  beam  has  been  given  by  Turpin6.  Non¬ 
coherent  implementation  using  an  intensity  modu¬ 
lated  light  emitting  diode  offers  inanunity  to  noise 
suffered  in  coherent  imaging.  The  light  budget  is 
an  important  design  consideration.  Source  collima- 
tion  is  required  for  efficient  acousto-optic  inter¬ 
action.  For  acousto-optic  modulation  at  low 
diffraction  efficiencies,  the  doubly  diffracted 
imaging  scheme  has  substantial  insertion  loss. 
Alternatively,  single  diffraction  imaging  may  be 
realized  by  utilizing  the  undiffracted  light  trans¬ 
mitted  by  A01  that  is  diffracted  into  the  first 
order  of  A02.  This  signal  is  mixed,  in  the  detec¬ 
tion  process,  with  the  first  order  signal  of  A01 
that  passes  undiffracted  by  A02.  Both  singly 
diffracted  signals  are  in  phase  and  may  be 
Schlieren  imaged.  Input  modulation  roust  be 
designed  such  that  the  undesired  terms  are  out 
of  band. 

The  next  examples  demonstrate  the  strength  of 
two-dimensional  detector  integration  processing. 

5 ,  AMBIGUITY  FUNCTION  PROCESSING 

Ambiguity  function  processing  is  important  in 
resolution  of  delay  and  Doppler  uncertainty.  In 
application  to  radar  signal  processing  or  external 
signal  parameter  measurement,  delay  and  Doppler 
are  often  time -varying.  An  ambiguity  function 
snapshot  is  therefore  desirable.  A  non-coherent 
optical  parallel  processing  implementation  for 
achieving  such  a  snapshot  is  given  in  this  section. 

This  example  utilizes  the  2-d  optical  system 
shown  in  Figure  5  which  was  generalized  in  the  last 
section.  Define  the  cross  ambiguity  function 
between  two  signals,  s^  and  s2,  as: 

X12<t,T)  '  4  Vt)S2(t'T)e~i2lIftdt  (U) 

where  T  and  f  are  the  delay  and  Doppler  variables. 
The  dual  input  optical  processor  output  is 
described  by  Equation  (12)  which  becomes: 


y<Vx2 
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B1 (t)e 


-iTTt 
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*  e  s1(t)s2(t-t1)e  *  dt 

(14) 

where  {a,b}  are  complex  chirp  waveforms.  This  is 
the  expression  for  the  cross  ambiguity  function 
with  a  quadratic  phase  term  that  may  be  cancelled 
through  post  detection  weighting.  The  real  archi¬ 
tecture  may  be  realized  by  the  method  described  for 
spectral  analysis,  in  fact,  the  ambiguity  function 
may  be  interpreted  either  as  a  multichannel  spec¬ 
trum  analyzer  given  by: 

X12(f.t)  -3\s1(t)s2(t-T)l  (15) 


for  multiple  time  delays,  I,  or  as  a  multichannel 
correlator  given  by: 

X12(£,t)  =  (Sl(t).e'l2lTft)*Is2(t))  (16) 


for  multiple  Dopplers,  f.  The  one-dimensional  pro¬ 
cessing  considerations  previously  discussed  are, 
therefore,  applicable.  The  Doppler  resolution  is 
commensurate  with  the  integration  time,  a  coarse 
resolution  may  be  maintained  during  acquisition 
period  and  a  higher  resolution  zoom  can  be  achieved 
by  longer  integration.  Multiple  correlation  peaks 
or  targets  can  be  processed  simultaneously  since  a 
linear  system  implementation  is  used. 

6.  TWO-DIMENSIONAL  SPECTRAL  ANALYSIS 

ft  non-coherent  optical  time  integrating 
approach  to  two-dimensional  complex  spectral  analy¬ 
sis  is  described.  This  method  is  an  extension  of 
the  technique  described  for  one-dimensional  analy¬ 
sis.  Very  large  bandwidth  (greater  than  10!) 
spectral  analysis  of  electrical  signals  can  be 
performed  using  this  time  domain  approach,  without 
requiring  storage  of  the  signal.  Spectral  analysis 
of  raster  scanned  video  imagery  can  be  performed  as 
well. 


The  optical  system  (Figure  5)  described  in  the 
last  sections  is  used  in  this  example  as  well.  The 
algorithm  is  implemented  with  three  reference  chirp 
waveforms.  The  source  is  modulated  with  the  input 
signal,  which  is  pre-multiplied  by  a  chirp,  and  the 
other  chirps  modulate  ftOl  and  A02.  The  chirp  rates 
are  chosen  such  that  the  product  of  the  three 
chirps  result  in  a  temporal  sine  wave,  with  fre¬ 
quency  that  is  varying  from  pixel-to-pixel  on  the 
detector  array.  Each  detector  element  integrates 
the  product  of  this  sine  wave  times  the  input 
signal,  and,  in  this  fashion,  produces  the  spec¬ 
trum.  The  frequency  difference  between  output 
lines  in  one  dimension  is  taken  to  be  N  times  the 
frequency  difference  between  lines  in  the  other 
dimension,  where  N  is  the  number  of  detectors  per 
line.  For  a  square  array,  the  number  of  spectral 
samples  is  proportional  to  N2 .  The  technology  is 
currently  limited  by  detector  array  size,  and 
detector  readout  rate. 

The  high  rate  chirps  are  folded  over  on  the 
order  of  N  times  during  the  integration  period, 
which  creates  a  comb  frequency  sampling  along  the 
coarse  frequency  axis,  with  comb  samples  being  a 


. ' 


‘function  of  the  fine  frequency  axis.  The  real 
implementation  is  performed  with  frequency  offset 
chirps,  similar  to  the  previous  examples.  Trans¬ 
forms  of  raster  scan  imagery  are  more  complicated, 
requiring  synchronization  that  is  not  described  in 
this  paper. 

High  resolution  spectral  analysis  by  detector 
integration  processing  has  several  advantages  over 
coherent  spatial  techniques.  These  include,  larger 
dynamic  range,  no  signal  storage,  complex  output, 
flexibility  in  processing,  advanced  device  tech¬ 
nology  (acousto-optic),  and  non-coherent  optical 
system  implementation. 

7.  SUMMARY 

Detector  integration  optical  signal  process¬ 
ing,  as  introduced  in  this  discussion,  has  been 
generalized  for  one-  and  two-dimensional  complex 
signal  processing.  These  techniques  are  imple¬ 
mented  by  time  integration  and  are  well  suited  to 
time  domain  algorithms.  The  optical  realization, 
utilizing  acousto-optic  input  devices  and  charge 
coupled  image  sensor  output,  creates  a  processor 
that  is  compatible  with  signal  processing  systems. 
Such  implementation  affords  very  large  time  band¬ 
width  signal  processing  without  the  storage  prob¬ 
lem  associated  with  coherent  spatial  techniques. 
Flexibility  of  the  transversal  filter  architecture 
is  achieved  by  means  of  variable  time  integration 
and  actively  generated  reference  waveforms.  The 
optical  system  described  is  multi-purpose  as 
illustrated  in  the  examples.  Complex  computation 
is  performed  without  requiring  a  coherent  optical 
system. 
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Abstract 

A  hybrid  optical  approach  to  coherent  spectrum  analysis  realizes  gains  associated  with 
time- integrating  techniques  while  still  retaining  some  of  the  advantages  of  space- 
integrating  techniques.  Measurement  of  large  time-bandwidth  spectra  is  Performed  with  high 
sensitivity,  efficiency,  and  range  of  detectable  signal  power  levels,  particularly  in  large 
bandwidth  applications,  where  an  optically  transformed  periodic  chirp  forms  a  distributed 
local  oscillator  as  a  reference  for  time  integration. 

Introduction 


The  continuing  advancement  of  electro-optical  component  technologies  is  increasingly  aid¬ 
ing  realization  of  optical  techniques  for  signal  processing.  These  developments  have  been 
encouraged  by  the  rapid  advancement  in  electronic  processing  fcechnioues  and  hardware.  This 
latter,  seemingly  paradoxical,  effect  derives  from  the  fact  that  the  progress  in  electronic 
processing  has  made  sophisticated  but  practical  processors  realizable  provided  the  burden 
of  very  high  throughput  rates  for  certain  operations  can  be  borne,  precisely  in  the  area  in 
which  the  parallel  processing  capabilities  of  optical  technioues  are  most  valuable. 

Spectrum  analysis  is  an  operation  that  has  become  particularly  well  addressed  bv  optical 
methods.  In  this  paper  we  will  focus  our  attention  on  two-dimensional  coherent  optical 
techniques  for  spectrum  analysis,  presenting  a  brief  review  of  the  snace-integratinn (D  and 
time- integrating '2' 3' 41  approaches  and  a  discussion  of  hybrid  implementations'5'  that  exhib¬ 
it  space- integrating  and  time-integrating  characteristics.  Our  attention  is  directed 
toward  the  examination  of  those  characteristics  that  make  optical  techniaues  most  useful  in 
processing  systems,  particularly  bandwidth  and  time-bandwidth  product,  and  a  comparison  of 
implementation  variables  such  as  optical  efficiency  and  signal  detectabilitv.  Since  soace- 
integrating,  time- integrating,  and  hvbrid  techniques  can  be  configured  in  two  dimensions  for 
compact,  large  time- bandwidth  processing,  the  discussions  below,  organized  to  emphasize  the 
close  parallel  among  the  techniques,  will  be  given  in  terms  of  two-dimensional  architectures. 

Coherent  Space-Integrating  and  Time-Integrating  Spectrum  Analysis 

Space- integrating  coherent  processors  rely  on  the  physics  of  coherent  light  propagation, 
which  permits  a  simple  lens  to  transform  a  light  amplitude  into  its  amplitude  spectrum. 

Since  this  can  be  performed  in  two  dimensions,  large  blocks  of  information  can  be  processed 
simultaneously.  Figure  1  shows  an  example'*'  of  a  two-dimensional  spectrum  analyzer  using 
acousto-optic  devices  to  deflect  and  modulate  the  laser  light.  Each  acousto-ootic  beam 
deflector  (AOBD)  causes  the  light  to  be  deflected  in  such  a  manner  that  the  combined  result 
is  that  the  light  focused  at  the  recording  medium  is  scanned  in  a  multiline  raster  as  it  is 
modulated  in  intensity  by  the  acousto-optic  modulator  (AOM) . 


Figure  1.  Two-dimensional  space- 
integrating  spectrum 
analyzer. 
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Since  the  recording  medium  responds  to  light  intensity,  the  signal  modulating  the  light 
power  must  be  on  a  bias.  The  recording  medium  is  used  in  a  linear  region  of  the  amplitude 
transmittance  versus  exposure  curve.  After  a  full  frame  has  been  recorded  L..e  medium  is 
illuminated  with  a  uniform  wave  of  coherent  light.  The  transmitted  amplitude  is  trans¬ 
formed  by  a  lens  to  produce  a  two-dimensional  intensity  distribution  that  is  proportional 
to  the  signal  power  spectrum  arranged  in  a  multiline  format,  one  coordinate  corresponding 
to  coarse  frequency  increments  and  the  other  to  fine  tuning;  in  other  words,  the  high  reso¬ 
lution  spectrum  is  folded  into  many  horizontal  lines  (in  Figure  1)  in  the  same  manner  that 
the  signal  is  folded  into  many  vertical  lines  on  the  recording  medium. 

An  analogous  two-dimensional  time-integrating  analyzer  (2)  is  shown  in  Figure  2,  config¬ 
ured  as  the  space-integrating  system  of  Figure  1,  except  that  the  recording  •  eaium  is 
removed  and  the  scanning,  modulated  light  is  transformed  by  a  lens  to  a  modulated  plane 
wave  deflecting  in  two  directions  at  a  detector  array  plane.  Here  it  interferes  with  a 
constant,  uniform  reference  wave  incident  from  an  oblique  angle.  The  signal,  without  a 


bias,  modulates  the  amplitude  of  the  light  at  the  AOM,  rather  than  the  intensity,  with  the 
result  that  the  time-integrated  energy  distribution  in  the  detection  plane  after  a  full 
frame  is  the  amplitude  spectrum.  This  technique  of  recording  the  spectrum  is  equivalent  to 
recording  an  off-axis  Fourier  transform  hologram  of  the  transparency  generated  in  the 
space-integrating  technique  of  Figure  1  while  scanning  it  with  a  focused  spot  rather  than 
illuminating  it  all  at  once,  as  in  Figure  1.  The  time-integrated  spectrum  is  an  amplitude 
spectrum  on  a  spatial  carrier,  still  maintaining  all  the  phase  information,  that  may  be 
electronically  demodulated.  The  exposure  in  the  detection  plane  of  the  space-integrating 
analyzer  is  a  power  spectrum,  requiring  no  demodulation;  if  an  amplitude  spectrum  is 
desired,  an  off-axis  reference  wave,  as  in  Figure  2,  may  be  added,  producing  a  Fourier 
transform  hologram  of  the  signal  transparency,  with  precisely  the  amplitude  spectrum  dis¬ 
tribution  as  that  of  Figure  2. 

The  time- integrating  technique  can  be  viewed  as  a  heterodyne  detection  technique  in 


DCTBCTION 

n*IH 


local  oscillator  refer¬ 
ence  waves. 
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which  the  signal  is  detected  by  comparison  with  a  local  oscillator  at  every  point  in  the 
detection  plane,  the  frequency  of  the  local  oscillator  being  a  function  of  the  coordinate 
of  the  point.  <2>  Figure  3  is  another  optical  scheme  for  a  time-integrating  analyzer;  while 
Figure  2  was  designed  to  illustrate  the  close  relationship  between  time-integrating  and 
space-integrating  techniques  (Figure  1),  the  rearrangement  of  Figure  3  is  meant  to  help 
clarify  the  concept  of  the  distributed  local  oscillator  by  creating  a  physical,  correspond¬ 
ing  to  the  mathematical,  separation  of  the  signal  and  the  local  oscillator  reference.  The 
relevant  term  in  the  power  density  integrated  over  time  T2  at  the  detector  olane  can  be 
written 


‘•yl 


S(t)  R (x,y, t)  dt. 


The  precise  form  of  the  local  oscillator  R(x,y,t)  is  discussed  below. 


Hybrid  Spectrum  Analysis 


The  above  two  approaches  to  high  resolution  spectral  analysis  exhibit  different  perfor¬ 
mance  bounds  and  implementation  requirements.  The  TI  technique  does  not  require  a  record¬ 
ing  material  for  intermediate  signal  storage  or  the  associated  dual  optical  system  (for 
recording  and  reading) .  Since  the  AOM  is  used  in  the  intensity  modulation  mode  in  the 
space-integrating  system,  it  has  only  half  the  bandwidth  capability  as  the  same  device  in 
the  time-integrating  system.  In  wide  bandwidth  applications  a  device  with  a  low  figure  of 
merit  (such  as  LiNb03>  roust  be  used  because  of  acoustic  attenuation.  However,  in  a  space- 
integrating  system  the  requirement  to  drive  the  AOM  with  a  bias  at  the  linear  region  of  the 
diffraction  efficiency  versus  drive  voltage  curve  may  be  impossible  to  meet.  On  the  other 
hand,  spatial  integration  is  capable  of  providing  much  higher  signal  levels  at  the  detec¬ 
tor,  and  the  power  is  localized  to  the  signal  area,  whereas  the  presence  of  a  signal  in  the 
time-integrating  scheme  contributes  a  bias  to  the  entire  detector.  The  significance  of 
this  is  that  the  presence  of  a  strong,  but  perhaps  uninteresting,  signal  does  not  interfere 
with  the  detection  of  a  much  weaker  signal  that  is  not  overlapping  in  frequency.  The  light 
power  incident  on  the  signal  transparency  can  be  increased  to  bring  a  weak  signal  into  the 
dynamic  range  of  the  detector,  allowing  the  strong  signal  to  saturate.  This  cannot  be  done 
with  a  time-integrating  system  because  a  strong  signal  adds  a  saturating  bias  to  the  entire 
detector  plane.  (The  important  considerations  of  dynamic  range  improvement  by  heterodyne 
detection  compared  with  direct  detection  do  not  qualify  this  observation  if,  for  comparison 
purposes,  a  constant  off-axis  reference  wave  is  added  to  the  space- integrated  spectrum  at 
the  detector. ) 


A  hybrid  spectrum  analysis  technique  that  succeeds  in  combining  the  above  important  fea¬ 
tures  of  space-integrating  and  time-integrating  architectures  is  illustrated  in  Figure  4. 
Two  important  features  distinguish  this  arrangment  from  that  of  Figure  3.  First,  the  opti¬ 
cal  aperture  of  the  AOM  has  been  increased  (to  become  an  AOBD,  by  our  definition) ,  thus 
generating,  by  spatial  integration,  an  amplitude  spectrum  at  the  detection  plane.  Second, 
the  reference  wave  has  been  redefined  to  be  the  Fourier  transform  of  the  wave  diffracted 
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from  the  AOBD  that  is  driven  by  the  fast  axis  periodic  chirp.  The  form  of  this  is  given 
below,  but  whatever  form  the  reference  wave  takes,  the  relevant  exposure  term  is 


E (x,y) 


A (x, t) 


R  ( x, y ,  t) 


dt  , 


(2) 


whe  re 


A  ( x ,  t )  =  I  e'l2,Tf't  LM  f '  )  sine  itT^avx-f')  df' 

J- oc 


(3) 


is  the  signal,  spatially  distributed  according  to  its  amplitude  spectrum,  the  estimate  of 
which  is  limited  by  the  delay  time  T^  of  the  AOBD.  v  is  the  acoustic  velocity,  a  =  2tt/\r, 
1  is  the  optical  wavelength  and  F  is  the  focal  length  of  the  transform  lens. 

Spatially  Distributed  Local  Oscillators 

The  reference  wave  of  Figure  4  can  be  described^*  by 


R (x,y , t) 


i  kx 


e 


£  e-i2irfp+(y/D)  ]  (t/Tj) 
P 


sine  irT  (A f  x/D)  -  (p/TjH , 


fx  <  p/T1  <  f  +Af 


(4) 


where  k  is  the  wave  vector  of  the  spatial  carrier,  f^  is  the  minimum  chirp  frequency  and  Af 
is  the  chirp  bandwidth.  (Unimportant  quadratic  ohase  factors  have  been  omitted  for  simpli¬ 
city.)  The  frequencies  f^  to  f^  +  Af  are  assumed  to  lie  in  the  Bragg  bandwidth  of  the  AOBD, 
and  the  chirp  period  is  equal  to  the  delay  time  Tj.  The  distributed  local  oscillator  des¬ 
cribed  by  Equation  4  and  illustrated  in  Figure  5  is  seen  to  consist  of  an  array  of  grid 
lines,  each  urid  line  having  width  1/T\  and  frequency  [p+ (y/D) ] /Tj . 


Fig-iv 


j . 


Representation  of  a 
two-dimensional  dis¬ 
tributed  local  oscil¬ 
lator  reference  ampli¬ 
tude.  A  periodic 
chirp  is  optically 
Fourier  transformed  ir. 
the  x  direction. 


Let  us  now  return  to  the  reference  wave  of  Figure  3  and  consider  the  effect  on  the  light 
amplitude  of  the  chirp  drive  for  the  fast  axis  only.  Considering  the  periodic  chirp  func¬ 
tion 


C:  (t) 


£  ei(a(t-nT^)2  +  b(t-nTj)] 
n 


rect  t-nT^, 

~fT~ 


a  =  irAf ,  b  =  2itf  i 
Tl 


(5) 


applied  as  a  Doppler  up-shift  by  the  AOM  and  a  down-shift  by  the  AOBD, 


the  net  result  is 


R  ( x ,  t ) 


£  eilbx  +  2axt-2anTi] 


rect  t-nTi  rect  t-nT^-x/v. 


(6) 


To  render  this  in  a  more  visible  form,  we  define  G(t-nTi)  and  use  the  identity 


G(t-nT,)  *  rect  t-nT,  rect  t-nTi-x/v  »  rect  x  rect  t-nTi-x/2v. 

A  —  ■  A  -  ■  •  “l.-fT 1  - ? - 


(7) 


2vT  i 


T^-  x/v 
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With  Equation  7  and  the  transform  relationship 


F(f) 


■/* 

•'-00 


i2nf  t 


G(t)dt  =  rect 


i2nf x/2v 


2vT, 


[tx-  |  v|]  sine  itf  [tj-IvI]  , 


(8) 


Equation  6  can  be  written,  after  some  manipulation,  as 

RU.t)  =  ^ibx  +  i2aXt/v£  .  JLje-i2*p$  -  ^]t  ( 

If  we  now  include  the  y  deflection,  neglecting  any  periodicity,  substitute  the  explicit 
form  for  F  (again  omitting  an  unimportant  quadratic  phase  factor) ,  and  include  an  angular 
offset,  the  two-dimensional  spatially  distributed  local  oscillator  becomes 


R(x,y,t) 


ikx  i2n[p+(y/D)]t/T1 


[i-IvT^l]  sinc{irTi[^’p][1-|^I,]}'0Sxsa  (10) 


Equation  10  represents  a  two-dimensional  distributed  local  oscillator  that  is  characterized 
by  degradation  in  high  frequency  resolution  when  the  fast  axis  AOBD  is  driven  at  high  chirp 
repetition  rates.  As  illustrated  in  Figure  6,  the  uniform  modulus  reference  wave  is  com¬ 
posed  of  a  set  of  grid  lines  (labeled  by  p) ,  with  spacing  l/Tj  having  frequencies 
lp+ (y/D) JT1(  as  does  the  distributed  local  oscillator  of  Eauation  4,  represented  in  Figure 
5.  However,  the  width  of  a  grid  line  is  larger  by  a  factor  [1- (x/VTi) J “1  and  the  peak  is 
lower  by  Il-tx/vTj^)].  Noting  that  the  optical  aperture  D  is  equal  to  |v|x,  where  t  is  the 
acoustic  transit  time  cf  the  cell,  the  line  width  degradation  factor  is  [1- (x/D) (x/Ti ) ] _1 . 
The  ratio  x/Tj  is  the  fraction  of  the  chirp  that  is  within  the  AOBD  aperture.  As  the  chirp 
period  T ^  approaches  x  the  grid  lines  get  broader.  This  spreading  results  in  an  ambiguity 
at  higher  frequencies  due  to  overlap  or  crosstalk  onto  neighboring  frequency  locations  in 
the  coarse  frequency  direction,  observed  as  ghost  signals  in  a  measured  signal  that  are 
separated  by  N2  resolvable  elements  (where  N2  the  time-bandwidth  product  in  the  y  direc¬ 
tion)  .  A  constraint  is,  therefore,  placed  on  the  chirp  period,  and  in  turn  on  the  band¬ 
width  and  resolution  of  the  system.  A  ratio  Tj/t  *  4  corresponds  to  a  33%  spread  at  the 
maximum  frequency,  or  about  10  dB  crosstalk  suppression  at  a  nearest  neighbor  peak.  If  we 
choose  this  to  be  the  maximum  allowable  error,  then  the  bandwidth  to  which  the  AOBD  must 
respond  (Nj/x)  is  four  times  larger  than  the  system  bandwith  (Nj/Tj) . 


Figure  6. 


Representation  of  a 
two-dimensional  dis¬ 
tributed  local  oscil¬ 
lator  reference  ampli¬ 
tude.  A  periodic 
chirp  is  imaged  in  the 
x  direction. 


System  Design  and  Performance  Considerations 

The  preceding  discussion  of  distributed  local  oscillators  serves  to  illustrate  one  of  the 
benefits  of  utilizing  the  Fourier  transform  Dronerties  of  lenses  with  soatially  modulated 
coherent  light.  Specifically,  the  full  bandwidth  capability  of  acousto-optic  devices  can  be 
exploited.  As  an  illustration  of  the  operational  limitations  involved  in  using  acousto¬ 
optic  cells,  Figure  7  contains  some  tradeoff  curves  relative  to  some  reoresentative  mater¬ 
ials.  The  solid  lines  represent  rather  soft  upper  bounds  to  the  time-bandwidth  product  as 
a  function  of  bandwidth,  assuming  a  device  bandwidth  to  center  freouenev  ratio  of  O.R,  and 
limited  by  acoustic  absorption.  An  acoustic  attenuation  of  3dB  over  the  optical  aperture 
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Figure  7.  Acousto-optic  device 
tradeoff  curves  for 
representative  mater¬ 
ials.  Solid  lines  de¬ 
fine  acoustic  attenua¬ 
tion  bounds.  (TeC>2)L 
and  (Te02)s  represent, 
respectively,  the  long¬ 
itudinal  and  slow  shear 
modes  of  Te02. 


At 


was  assumed  here.  The  dashed  lines  are  lines  of  constant  cell  length;  30mm  reoresents  a 
practical  upper  limit  in  most  cases.  Figure  7  does  not  nretend  to  summarize  all  the  con¬ 
straints  involved  in  fabricating  acousto-ootic  devices,  even  for  the  materials  reDresented, 
but  to  provide  a  reasonable  basis  of  discussion  of  the  operating  regions  of  devices  that 
might  be  realized  for  the  systems  we  are  presently  considering. 

As  an  example,  a  LiNbOg  AOBD  with  1  GHz  bandwidth  and  time-bandwidth  product  of  103  is 
within  the  realm  of  existing  technology,  and  may  be  used  to  implement  a  1  GHz  spectrum  ana¬ 
lyzer  only  if  the  distributed  local  oscillator  is  of  the  spatially  transformed  chirp  varie¬ 
ty  in  the  fast  axis  direction  (Equation  4) .  The  number  of  resolvable  frequencies  of  such  a 
system  would  be  lo6  if  the  slow  axis  factor  to  the  two-dimensional  local  oscillator  was 
implemented  as  an  imaged  AOBD  (rather  than  optically  Fourier  transformed) ,  using  a  25mm  Te02 
device  in  the  slow  shear  mode  with  a  chirp  bandwidth  of  25  MHz.  The  slow  chirp  period  would 
be  1  ms,  so  that  T2/t  =  25  for  the  slow  axis  device,  large  enough  to  avoid  significant  chirp 
fill  time  effects. 


As  a  general  rule  the  diffraction  efficiency  of  an  AO  device,  for  the  maximum  rated  drive 
power,  decreases  with  the  bandwidth  of  the  device,  at  least  in  the  high  frequency  region. 
Requiring  a  device  several  times  wider  in  bandwidth  than  the  signal  to  be  analyzed  could  be 
a  significant  factor  in  optical  power  considerations.  This  is  especially  true  where  a  change 
in  acousto-optic  material  is  involved.  For  example,  a  120  MHz  system  with  700  resolvable 
frequencies  in  the  coarse  direction  could  be  implemented  with  a  25mm  Te02  device  operated  in 
the  longitudinal  mode  in  a  hybrid  system.  A  purely  time-integrating  system  would  require 
about  500  MHz  of  device  bandwidth,  imolving  a  LiNbOg  AOBD,  a  material  with  a  much  lower  fig¬ 
ure  of  merit,'61  and,  as  a  result,  lower  optical  efficiency. 

Hybrid  systems  retain  an  important  characteristic  not  shared  with  ourely  time-integration 
spectrum  analysis;  very  weak  signals  can  be  detected  in  the  presence  of  strong  signals  not 
overlapping  in  coarse  frequency  lines  by  raising  the  laser  power  incident  on  the  signal 
AOBD,  permitting  the  st  ong  signal  to  saturate  its  line  and  raising  the  power  of  the  weak 
signal  into  the  dynamic  range  of  the  detector.  In  addition  to  the  optical  efficiency  gains 
resulting  from  device  bandwidth  considerations,  a  large  gain  in  efficiency  in  the  signal 
path  of  the  hybrid  system  of  Figure  4  is  achieved  by  concentrating  the  signal  light  into  the 
spectrum,  rather  than  distributing  it  uniformly  over  the  detector. 

The  coherent  hybrid  approach,  like  the  purely  time-integrating  technique,  avoids  the 
necessity  of  using  an  intermediate  recording  medium  and  associated  optical  readout  com pli¬ 
cations.  On  the  other  hand,  such  recording  media  could  be  used  in  the  detection  Diane  to 
produce  an  intermediate  recording  of  the  amplitude  spectrum,  with  real-time  optical  demodu¬ 
lation.  This  optical  demodulation  by  coherent  illumination  and  spatial  filtering  provides 
the  two-dimensional  detector  array  with  the  power  spectrum,  decreasing  the  number  of 
required  array  elements  becaus^.of  the  absence  of  the  carrier,  and  permitting  multi-frame 
time-integration  detection.'  '  (It  should  be  emphasized  that  when  recording  media  are 
used  in  the  detection  plane  of  hybrid  systems  the  bandwidth  advantages  of  such  architectures 
are  retained.)  There  is  a  degree  of  versatility  of  detection  techniaues  therefore  associa¬ 
ted  with  time-integrating  and  hybrid  schemes,  a  fact  that  may  be  of  some  significance  in  the 
wide  bandwidth  applications  we  are  addressing,  where  the  state  of  the  art  in  detection  hard¬ 
ware  is  presently  the  limiting  implementation  problem. 
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Summary 

The  hybrid  approach  to  the  two-dimensional  coherent  spectrum  analysis  combines  some 
important  features  of  space-integrating  and  time-integrating  techniques.  Spatially  trans¬ 
formed  periodic  chirps  provide  a  distributed  local  oscillator  that  permits  time-integration 
detection  at  the  full  bandwidth  of  acousto-optic  devices.  The  range  of  detectable  signal 
levels  can  be  extended  well  beyond  the  instantaneous  dynamic  range  of  the  detector  because 
of  the  spatial  separation  of  the  light  energy  in  the  spectrum.  This  spatial  spectrum  also 
results  in  higher  optical  efficiency;  overall  efficiency  is  improved  further  by  permitting 
the  use,  for  the  same  system  bandwidth,  of  acousto-optic  devices  of  lower  bandwidth  and 
therefore,  higher  achievable  diffraction  efficiency,  a  reflection  of  the  tradeoffs  involved 
in  the  performance  of  such  devices.  Some  versatility  exists  in  the  methods  available  for 
optical  detection,  the  aspect  that  presents  the  most  difficult  implementation  problem  for 
wide-bandwidth  systems. 
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Time  integrating  optical  signal  processing 

P.  Kellman  Abstract.  Time  integrating  acousto-optic  processors  realize 

flexible,  multipurpose  complex  signal  processing  architectures 
ESL,  Incorporated  based  on  correlation  algorithms.  One-  and  two-dimensional 

495  Java  Drive  techniques  are  presented  including  examples  of  spectral 

Sunnyvale,  California  94086  analysis  and  ambiguity  function  processing.  Noncoherent  opti¬ 

cal  processor  implementation  using  interferometric  detection 
with  electronic  reference  is  described  and  experimental  results 
are  given. 


Keywords,  acousto-optics .  signal  processing,  correlators,  spectra!  analysis,  ambiguity 
function,  triple  product  processor,  chirp  algorithm. 
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I.  INTRODUCTION 

Optical  techniques  tor  linear  signal  processing  lend  themselves 
naturally  to  large  time-bandwidth  product  operations  due  to  the 
high  degree  of  parallelism  found  in  optical  systems.  Parallel  optical 
signal  processing  has  traditionally  exploited  spatial  integration1  1 
to  realize  tunctions  such  as  filtering  and  spectral  analysis.  In  this 
manner,  the  potential  time-bandwidth  product,  or  number  of 
cycles  integrated,  is  proportional  to  the  large  number  of  degrees  of 
freedom  ot  the  optical  system,  though  a  practical  limitation  is 
imposed  by  input  and  output  devices  such  as  light  modulators  and 
detectors. 

Integration  in  time  rather  than  space  may  be  used  to  realize  ex¬ 
tremely  large  time-bandwidth  products.  Time  inte^  ating  techni¬ 
ques  have  recently  been  generalized  for  one-  and  two-dimensional 
complex  signal  processing.'"1  Interest  in  time  integrating  techniques 
has  resulted  due  to  the  attractive  device  technology  as  well  as  flexi¬ 
bility  of  time  integrating  algorithms.  An  important  consequence  of 
time  integrating  techniques  is  the  ability  to  operate  on  signals  with 
very  large  time-bandwidth  product  without  having  to  store  the 
entire  time  history  as  a  spatial  record.  Therefore,  an  important 
class  of  two-dimensional  and  multichannel  processing  algorithms 
may  be  performed  without  requiring  two-dimensional  spatial  light 
modulators.  Acousto-optic  devices  and  charge  coupled  image  sen¬ 
sors  are  particularly  well  suited  for  time  integrating  processor 
implementation.  Further,  these  signal  processing  techniques  may  be 
realized  with  either  coherent  or  noncoherent  optical  systems. 

The  optical  signal  processing  architectures  discussed  in  this 
paper  are  based  on  acousto-optic  or  other  traveling  wave  input 
modulation  devices.  Time  integrating  optical  correlators  were  first 
demonstrated  using  translating  optical  masks,*"’  and  scanning 
detector  systems.’  10  In  signal  recording  applications,  time  inte¬ 
grating  techniques  have  been  used  to  compensate  Doppler  shift.11  11 
Time  integrating  correlator  implementation  has  been  demonstrated 
using  acousto-electric  surface  wave  technology.1'  Acousto-optic  im¬ 
plementations  of  one-dimensional  time  integrating  correlation  and 
spectral  analysis  have  first  been  introduced  by  Montgomery,1* 
Sprague,  and  Koliopoulos.1’  Two-dimensional  time  integrating 
techniques  were  introduced  by  Kellman’  *  and  Turpin.'  In  this 
paper,  time  integrating  acousto-optic  signal  processing  is  reviewed, 
and  the  concept  of  interferometric  detection  with  electronic  reference 
is  described.  These  techniques  are  generalized  for  complex  computa¬ 
tion.  and  examples  of  spectral  analysis  and  ambiguity  function  pro¬ 
cessing  are  given. 

This  paper  is  organized  as  follows.  Time  integrating  correlation 
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is  reviewed  in  Section  2  and  interferometric  implementation  with 
electronic  reference  is  described.  Time  integrating  spectral  analysis 
by  means  of  the  chirp  algorithm  is  described  in  Section  3.  and  the 
experimental  result  of  a  noncoherent  optical  implementation  is 
given.  Two-dimensional  processing  is  presented  in  Section  4  and 
examples  of  ambiguity  function  processing  and  spectral  analysis  are 
given  in  Sections  5  and  6. 

II.  TIME  INTEGRATING  CORRELATION 


The  one-dimensional  time  integrating  correlator  is  reviewed.  Con¬ 
sider  the  optical  realization  shown  conceptually  in  Figure  I .  The 


light  source  is  temporally  modulated  by  a  signal,  s,(t),  to  produce 
an  output  intensity  lift),  which  illuminates  an  acoustic  delay  line 
light  modulator.  A  second  signal,  silt),  is  introduced  into  the 
acoustic  cell,  which  spatially  modulates  the  source  intensity  by  an 
amount  denoted  by  I2(t  —  x/v),  where  x  is  the  spatial  dimension  and 
v  is  the  acoustic  velocity.  The  acoustic  signal  is  schlieren  imaged 
onto  a  linear  photodiode  array;  the  intensity  distribution  in  the 
image  plane  is  given  by  the  product  Ii(t)l2(t  —  x/v).  Distinction  be¬ 
tween  coherent  and  noncoherent  optical  system  implementation 
and  interferometric  versus  noninterferometric  detection  will  be 
described. 

The  charge  integration,  directly  proportional  to  exposure,  is 
performed  by  detectors  at  discrete  positions,  x;.  The  resultant  out¬ 
put  voltage  at  the  i'*1  detector  element  is  therefore; 

R|l<Tj)  =  f  Ii(t)Ij(t  — Xj/V)dt  (1) 

JT 

where  rj  —  xj/v,  and  T,  the  integration  time  of  the  detector,  is  set 
by  the  timing  of  the  charge  transfer  readout  register.  This  is  the 
desired  correlation,  where  I,  and  12  are  directly  related  to  the  input 
signals  s,  and  s2.  Several  modulation/detection  schemes  may  be 
employed  in  order  to  achieve  a  linear  relationship. 

The  range  of  relative  delay  between  input  signals  is  limited  by 
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the  acoustic  delay  length,  and  the  integration  time  is  set  by  the 
detector  readout  period  The  correlator  time-bandwidth  product, 
BT.  may  therefore  be  larger  than  the  delay  line  time-bandwidth 
product,  Bt,  since  the  integration  time  period  may  be  longer  than 
the  acoustic  delay.  In  spatial  integrating  correlators, '*  the  integra¬ 
tion  period  is  determined  by  the  acoustic  delay;  in  configurations 
employing  fixed  reference  masks,  the  range  of  relative  delay  is 
unlimited,  however,  in  correlators  with  acoustic  reference,  the 
range  is  limited  to  the  acoustic  delay.  Long  integration  and  flexibil¬ 
ity  of  variable  integration  are  achieved  with  the  time  integrating 
approach. 

The  achievable  processing  gain  is  limited  by  the  correlator 
dynamic  range  which  is  determined  by  the  detector  dynamic  range 
and  the  signal-to-bias  ratio.  The  signal-to-bias  ratio  depends  on  the 
light  intensity  modulation  depth.  Detector  dynamic  range  is  limited 
by  saturation  and  is  defined  as  the  ratio  of  saturation  level  to  rms 
noise.  The  processing  gain  and  dynamic  range  can  be  extended, 
however,  by  post-detection  digital  integration. 

In  all  approaches  to  time  integrating  optical  correlation,  the  goal 
is  to  achieve  a  term  in  the  detected  light  intensity  that  is  propor¬ 
tional  to  the  product  of  the  input  signals.  Acousto-optic  devices 
modulate  optical  phase,  thus  are  basically  nonlinear  modulators  of 
electric  field  amplitude  or  intensity.  However,  linear  electric  field 
modulation  is  approximated  at  low  diffraction  efficiency  (depth  of 
phase  modulation).  Linear  intensity  modulation  is  approximated 
by  operation  with  an  optical  bias  at  high  diffraction  efficiency  (*74 
phase  shift)  and  a  small  signal  modulation  depth.  This  latter  techni¬ 
que  has  been  exploited  by  Sprague  and  Koliopoulous.11  Inter¬ 
ferometric  detection  may  be  used  when  acousto-optic  modulation 
of  electric  field  is  linear.  In  one  implementation,  a  coherent  optical 
reference  beam  is  used  for  interferometric  detection.  Another  im¬ 
plementation  is  described  that  uses  an  electronic  reference  to  realize 
interferometric  detection  with  either  coherent  or  noncoherent  light. 
This  approach  is  realized  by  adding  a  reference  oscillator  signal  to 
the  acoustic  modulation,  rather  than  using  a  reference  beam  with 
separate  optical  path  The  basic  difference  between  the  coherent 
and  noncoherent  implementation  is  the  spatial  and  temporal  diver¬ 
sity  of  the  illumination.  Tolerance  to  angular  and  wavelength 
dispersion  is  determined  by  the  acoustic  diffraction  and  the  Bragg 
condition.  The  magnitude  transfer  function,  MTF,  is  related  to  the 
illumination  and  acoustic  field  by  convolution  of  their  angular 
spectrums. 

In  the  interferometric  schemes,  it  is  assumed  that  acousto-optic 
modulation  of  the  electric  field  is  linearly  proportional  to  the  drive 
voltage,  and  that  the  imaging  optics  pass  only  the  first  diffraction 
order.  This  condition  is  approximated  at  low  diffraction  efficiency. 
Third  order  intermodulation  may  be  in  band  and  may  contribute  to 
the  output.  In  special  cases,  these  terms  time-integrate  to  zero. 

The  implementation  described  uses  an  internally  modulated 
diode  source  as  shown  in  Figure  1.  A  reference  oscillator  signal  is 
added  to  the  acousto-optic  deflector  input  with  a  frequency  that  is 
offset  from  the  signal  modulation;  in  addition  the  illumination 
source  must  be  modulated  on  a  carrier  with  equal  frequency  offset. 
Both  single  and  double  sideband  modulation  are  analyzed. 

The  image  intensity  distribution,  1,  is  given  by  the  product  of 
the  source  modulation,  I,,  and  acousto-opftic  modulation  denoted 
by  I,  -  Ej  ‘  where  E,  is  the  complex  electric  field  modulation. 

lit  1 1  -  I,(t)I,(l  -  x/ v)  .  (2) 

for  double  sideband  modulation 

I  n>  -  A , 1 1  v2miS,(ti  cos|2*fot)J 

f  .  t  -  A)  *|l  ♦  v2m,s,(t)e je^^c' 

I. <t  -  F,'t>  ' 

A, |  J  *  *■  2\ftm,s,<t>cos|2»ft,l)) ,  (3) 

frequency  difference  between  the  reference 


oscillator  at  fc  and  the  double  sideband  suppressed  carrier  modula¬ 
tion  at  fc  +  f„.  The  +  1  diffraction  order  is  passed  by  the  imaging 
optics.  A,  and  A,  correspond  to  light  intensity  and  diffraction  effi¬ 
ciency  respectively;  m,  and  mj  are  constants  that  determine  the 
modulation  depth.  It  is  assumed  that  signals  s,  and  s,  are 
bandlimited  to  a  bandwidth  B  (i.e.,  |S(f)|  «  0,  |f|>  B).  and  have 
unit  average  power.  The  spectrum  of  the  input  modulation  is 
shown  in  Figure  2.  The  device  bandwidth  is  fQ  +  B. 


(c)  acousto-optic  deflector  modulation  for  DSB  example. 

For  fQ  >  3B  several  cross  terms  effectively  integrate  to  zero  and 
the  output  becomes  approximately 

R(t)=  AiA2[T  +  2m|*  jsj(t  —  r)dt 

+  2mlm,cos[2xf0T]*  js,(t)s,(t-r)dtj  (4) 

where  r  *  x/v.  The  first  two  terms  are  bias  and  the  last  term  is  the 
desired  correlation  on  a  spatial  carrier,  fc.  The  bias  terms  may  be 
eliminated  through  filtering.  The  ratio  of  signal-to-bias  for  maxi¬ 
mum  correlation  is  given  by  0  *  2m,mj/(l  +■  2m{), 

For  single  sideband  (SSB)  modulation  the  bandwidth  require¬ 
ment  is  cut  in  half.  The  spectrum  of  the  input  modulation  is  shown 
in  Figure  3.  For  fD  >  3/2  B  the  output  correlation  is  approximately: 

R(r)  =  A,Ai^jdt  +  m5j[sJ(t- r)  +  rj(t  — r)]dt 
+  2m,mi(Ru(T)cos!2*foT]+  R  u(T)sin|2Tf0rl)j  (5) 


Figure  3.  Spectrum  of  (a)  diode  source  modulation,  and  (b)  acousto¬ 
optic  deflector  modulation  for  SSB  example. 


where  Rjj(t)  is  the  desired  correlation  and  R  u  is  the  Hilbert 
transform  of  R„(R  u  “  Rifb  The  desired  correlation  R,,  may  be 
synchronously  detected. 

Complex  correlation  using  real  computation  is  described  in  the 
following.  In  general,  calculation  of  complex  multiplication  or  cor¬ 
relation  requires  four  real  operations.  Alternatively,  frequency 
translation  may  be  utilized  to  realize  complex  correlation  with  a 
single  real  correlator  at  the  expense  of  bandwidth.  The  time  inte- 
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grating  spectrum  analyzer  described  in  the  next  section  is  an  exam¬ 
ple  o(  complex  correlation  by  this  method.  Further  advantage  of 
correlation  at  a  carrier  frequency  is  the  ability  to  filter  the  desired 
correlation  from  additional  bias  terms.  The  expressions  for  com¬ 
plex  correlation  are  derived  in  the  following. 

Define  the  correlation  between  complex  signals  x(t)  and  y(t)  by: 

Rxy(r)  =  <  x(t)y  *(t  —  t)  >  (6) 

where 

x(t)  =  xR(t)  +  ix|(t) 

y(t)  =  yR(t)  +  iy,(t) 

and  <  •  >  denotes  ensemble  average. 

Re|Rxy<T)|  =  <xR(t)yR(t-r)>  +  <x,(t)y,(t  -  t)  > 

lm|Rxy(r)|*=  —  <xR(t)yI(t-r)>  +  <X|(t)yR(t-r)>  .  >7) 

Consider  the  real  correlation  between  x0(t)  and  y0(t). 

RXoyo  (r)  =  < x0(t)y0(t  —  x) > 

x0(t)  =  Re|x(t)e|2lf“‘  |  (8) 

Volt)  *  Relydle12**”1  | 

where  x(t)  and  y(t)  are  bandlimited  with  bandwidth  B,  and  f„  >  B. 

Rx,y„(^)  =  T-RelRxy(T)eiW1-  <9) 

Cross  products  at  the  sum  frequency  2f0  average  to  zero  for  f0  > 
B.  Thus,  the  correlation  between  x0  and  y„  is  at  a  carrier  frequency 
f0.  with  carrier  phase  modulated  by  the  phase  of  RXy,  and  envelope 
equal  to  the  magnitude  of  RXy.  The  real  and  imaginary  parts  of 
RXy(r)  may  be  derived  from  RXoy0  by  synchronous  detection  (in 
quadrature). 

III.  TIME  INTEGRATING  SPECTRAL  ANALYSIS 

An  approach  to  real  time  optical  spectral  analysis  of  electrical 
signals  that  uses  a  time  integrating,  rather  than  spatial  integrating, 
version  of  the  chirp  z-transform.  is  described.  The  large  time- 
bandwidth  correlation  is  accomplished  by  means  of  the  correlator 
described  in  the  previous  section.  A  time  integrating  architecture 
for  correlation  and  spectral  analysis,  using  oppositely  traveling 
acoustic  waves,  has  been  described  by  Montgomery.1*  For  time 
integrating  realizations,  frequency  resolution  is  determined  by  the 
detector  integration  period,  rather  than  the  acoustic  delay  length, 
therefore,  higher  resolution  can  be  achieved  than  by  spatial  integra¬ 
tion.  Variable  resolution  is  attainable  through  variable  time  inte¬ 
gration.  Use  of  intensity  modulation  leads  to  detection  of  the 
magnitude  spectrum  rather  than  the  power  spectrum  which  yields  a 
considerable  increase  in  dynamic  range.  This  approach  to  spectral 
analysis  achieves  a  flexibility  not  readily  achieved  by  coherent  opti¬ 
cal  spatial  methods.  The  chirp  algorithm,1’  "  is  reviewed,  and  the 
real  implementation  is  discussed. 

The  Fourier  transform  integral 

S(f)  =  jslUe*'2*'  dt  (10) 

may  be  rewritten  as 

S(f)  -  e-'**’ js(t)e"'I,‘elr<,'f)1  dt  .  (11) 

by  expressing  f«t  as  -Ft1-  (t  -f)'|  The  complex  realization  is 
shown  in  Figure  4.  The  post-phase  weighting  may  be  ignored  if 
only  the  magnitude  or  power  spectrum  is  required.  This  realization 


Figure  4.  Complex  realization  of  chirp  algorithm. 


converts  the  integral  transform  to  time  invariant  filtering  of  a 
preweighted  signal,  followed  by  post-weighting. 

The  chirp  algorithm  may  be  realized  with  a  single  real  correla¬ 
tion  by  translating  the  real  chirp  to  a  carrier  frequency,  f0,  greater 
than  the  signal  analysis  bandwidth,  B.  The  desired  difference  fre¬ 
quency  is  translated  to  baseband,  and  the  sum  frequency  is  an  out- 
of-band  chirp  that  integrates  to  zero.  The  product  of  the  difference 
frequency,  which  is  a  linear  function  of  delay,  and  the  input  signal, 
is  integrated,  producing  the  Fourier  transform.  The  magnitude 
spectrum  corresponds  to  the  detected  envelope.  The  phase 
modulates  the  spatial  carrier  and  may  be  detected  as  well,  resulting 
in  a  complex  output. 

Frequency  resolution  of  this  spectrum  analyzer  is  inversely  pro¬ 
portional  to  the  detector  integration  time,  T.  The  resolution  is  not 
set  by  the  acoustic  time  delay  which  typically  limits  the  use  of 
acousto-optical  processing  to  wideband  analysis.  A  wide  range  of 
integration  times  are  easily  realized  (0.1-100  msec)  using  self- 
scanned  arrays  of  moderate  size  (e.g.,  1024).  Integration  time 
periods  can  be  extended  by  means  of  electronic  accumulation. 

The  delay  line  length,  t,  and  the  number  of  detectors,  N,  set  the 
delay,  A r,  between  detector  elements  in  the  image  plane.  This 
determines  the  chirp  bandwidth  B'.  The  required  resolution  or  inte¬ 
gration  period  fixes  the  chirp  rate,  a,  (see  Figure  5).  The  spatial 


bandwidth  equals  the  chirp  bandwidth,  B',  plus  carrier  frequency, 
f0.  Let  k  be  the  number  of  samples  (detectors)  per  cycle  of  the 
highest  frequency  (k  —  2  for  Nyquist  sampling),  and  define  k'  as 
the  ratio  of  frequency  offset  to  chirp  bandwidth,  f0/B'.  The  above 
relationships  may  be  written  as: 

N-  k(f„  +  B')r 
a  -  B7T  “  B/r 

k'  -  f„/B’  £  1  .  (12) 

The  time-bandwidth  product  is,  therefore 

BT  -  N/k(l  +  k')  <  N/4  .  (13) 
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The  dynamic  range  is  limited  by  detection  noise  It  there  are  M 
distinguishable  levels  in  the  output  voltage  between  the  noise  floor 
and  the  photo  element  saturation  value,  and  if  0  is  the  maximum 
ratio  of  signal-to-bias,  then  the  dynamic  range  in  dB  is  given  by  20 
logm  M0/(1  +  0),  since  incident  exposure  is  directly  proportional  to 
spectral  magnitude.  This  results  in  a  considerable  increase  in 
dynamic  range  as  compared  to  coherent  optical  power  spectral 
analysis  in  which  incident  exposure  is  related  to  spectral  power 
density  and  the  dynamic  range  is  given  by  10  log,0M. 

An  experimental  breadboard  was  constructed  using  a  non¬ 
coherent  implementation  with  a  light  emitting  diode  source.  The 
modulation  was  double  sideband  and  the  chirps  had  octave  band¬ 
width  (k'  =  l).  The  correlator  output  fora  sine  wave  input  is  shown 
in  Figure  6.  In  this  photograph  the  envelope  is  observed  since  the 


Figure  6.  Time  integrating  spectrum  analyzer  output. 


exposure  covers  multiple  integration  periods  in  which  the  signal 
phase  is  changing.  Parameters  of  this  example  were:  N  =  1000, 
t"=40  gsec.  B=2.5  MHz,  T  =  2  msec,  BT  =  100,  k  =  5,  k'  =  l.  The 
signal-to-bias  ratio,  0.  was  20%  and  the  detector  dynamic  range 
(peak  to  rms )  was  M  =  1000.  The  dynamic  range  observed  was  40 
dB  (UM/2)  since  the  bias  was  set  at  M/2  rather  than  optimized  at 
M/(I+£).  The  dynamic  range  is  3  dB  greater  using  SSB  input 
modulation. 

IV.  TWO-DIMENSIONAL  PROCESSING 

Two-dimensional  integral  transforms,  for  which  the  kernel  is 
decomposable  in  the  proper  way,  may  be  implemented  by  time 
integrating  optical  processing.  Two  important  examples  are  am¬ 
biguity  plane  processing  and  spectral  analysis. 

Consider  the  optical  realization  shown  conceptually  in  Figure  7. 
The  optical  train  has  a  modulated  illumination  source,  two  acoustic 
delay  line  light  modulators,  and  a  matrix  array  of  detectors.  This 
configuration  may  be  employed  for  several  functions,  determined 
by  input  signal  and  reference  waveforms.  It  is  assumed  throughout 
this  discussion  that  two-dimensional  processing  is  applied  to  either 
very  long  one-dimensional  signals  or  to  two-dimensional  signals 
(e  g.,  imagery)  that  are  in  raster  format.  In  this  way,  the  input  will 
be  a  function  of  one  variable,  t;  the  output  R(r,,rj)  is  a  function  of 
two  variables. 

The  light  diffracted  by  the  first  acousto-optic  modulator  A01,  is 
diffracted  by  the  second,  A02,  in  an  orthogonal  direction.  Both 
acoustic  signals  are  imaged  on  the  detector  plane.  The  desired 
image  plane  intensity  distribution  has  a  term  proportional  to  the 
product  of  acoustic  signals,  s,(t  —  r, )s,(t  —  Tj),  where  r,  — y/v  and 
r,  —  x/v  are  time  variables.  The  resultant  detected  output  voltage  is 
proportional  to  the  integrated  charge  as  in  the  previously  described 


Figure  7.  Two-dimensional  time  integrating  optical  processor. 


correlator.  This  structure  implements  integrals  of  the  form: 

R(7i,T,)  =  |  Soft )S)(t  Ti )S2(t  7j)dt  .  (14) 

In  this  configuration,  all  signals  are  real.  In  order  to  perform 
complex  computation,  similar  considerations  apply  that  were  given 
for  one-dimensional  processing. 

Hybrid  realizations  that  use  a  combination  of  spatial  and  time 
integration  can  be  used  to  implement  a  variety  of  other  operations. 
An  example  of  spectral  analysis  utilizing  spatial  integration  for 
coarse  resolution  and  time  integration  for  fine  resolution  is  des¬ 
cribed  by  Bader.” 

The  two-dimensional  time  integrating  optical  processor  may  be 
implemented  with  several  modulation/detection  schemes.  A 
description  of  a  coherent  optical  implementation  using  an  inter¬ 
ferometric  optical  reference  has  been  given  by  Turpin.'  Non¬ 
coherent  implementation  using  an  intensity  modulated  light  emit¬ 
ting  diode  offers  immunity  to  noise  suffered  in  coherent  imaging. 
Noncoherent  optical  implementation  using  a  reference  oscillator 
for  interferometric  detection  may  be  used  to  realize  complex  opera¬ 
tion.  Input  modulation  must  be  designed  such  that  undesired  terms 
are  out  of  band. 

The  next  examples  demonstrate  the  strength  of  two-dimensional 
time  integrating  processors. 

V.  AMBIGUITY  FUNCTION  PROCESSING 

Ambiguity  function  processing  is  important  in  resolution  of  delay 
and  Doppler  uncertainty.  In  application  to  radar  signal  processing 
or  external  signal  parameter  measurement,  delay  and  Doppler  are 
often  time-varying.  An  ambiguity  function  snapshot  is  therefore 
desirable.  A  noncoherent  optical  parallel  processing  implementa¬ 
tion  for  achieving  such  a  snapshot  is  given  in  this  section. 

This  example  utilizes  the  2-D  optical  system  shown  in  Figure  7 
which  was  generalized  in  the  last  section.  The  cross  ambiguity  func¬ 
tion  between  complex  signals  x(t)  and  y(t)  may  be  defined  as 

{T 

x(t)y*(t  —  dt  (15) 

0 

where  the  variables  r  and  f  may  be  interpreted  as  delay  and  Dop¬ 
pler,  respectively.  A  real  implementation  of  the  complex  cross¬ 
ambiguity  function  is  described  that  uses  the  architecture  of  Section 
4.  The  spectral  analysis  operation  is  performed  using  the  chirp 
algorithm.  Complex  correlation  was  derived  in  Section  2  and  is 
easily  extended  to  two-dimensional  computation. 

Define  the  complex  input  signals  s„(t),  s,(t),  and  s2(t)  by 

s„(t)  -  x(t)e‘i<,wt' 
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Silt)  =  y*(t) 

s,(t)  -  emw,‘  ,  (16) 

and  consider  the  interferometric  implementations  of  Section  2  with 
electronic  reference.  In  the  case  of  complex  inputs,  the  intensity 
modulation  is  given  by 

I„(t)  =■=  A„[]  +•  v/2m<,Rejs0(t)e‘2’<“,  j  ] 

I,(t)  =  A,[l  +  2mjls,(t)!'  +  2v^miRe|s,(t)et2**,t  j  ] 

lift)  =  A2[i  +  2rrij  |s2(t)  ( 2  +  2v/2mJRe{s,(t)e"i2l<I,J  ]  (17) 

where  double  sideband  modulation  is  assumed.  The  integrated 
intensity.  R(r,,Tj), 

Rtr.-r,)  =  (  l0(t)l,<t-r,)Mt-r,)dt  (18) 

contains  the  desired  ambiguity  function  term. 


DELAY  AOuitc 


Figure  8.  Ambiguity  function  for  short  code. 


Rej  Axy(T,,uT2)e~'l2l(,,r'  "  f*,‘)  +  “"i1)  ,  (19) 

which  is  at  a  carrier  frequency  and  may  be  filtered  from  other 
cross-product  terms.  The  quadratic  phase  term  may  be  cancelled 
through  post-detection  weighting.  The  variables  r,  and  r,  corres¬ 
pond  to  delay  and  Doppler  respectively.  The  Doppler  resolution  is 
commensurate  with  the  integration  time;  the  analysis  bandwidth  is 
determined  by  Eq.  (13).  A  coarse  resolution  may  be  maintained 
during  a  signal  acquisition  period,  and  a  higher  resolution  room 
can  be  achieved  by  longer  integration.  Multiple  correlation  peaks 
or  targets  can  be  processed  simultaneously  since  a  linear  system 
implementation  is  used. 

Ambiguity  function  processing  was  demonstrated  using  the 
noncoherent  optical  implementation  shown  in  Figure  7.  Devices  in¬ 
cluded  a  Hitachi  HLP-20  light  emitting  diode,  Fairchild  SL62926 
charge  coupled  device  image  sensor,  and  Isomet  acousto-optic 
devices.  The  diode  has  a  30  MHz  3  dB  bandwidth  and  was  biased  at 
an  average  optical  power  of  approximately  10  mW.  The  image  sen¬ 
sor  has  380  X  488  elements  and  was  operated  with  an  integration 
period  of  33.34  msec.  The  acousto-optic  devices  have  a  30  MHz  1 
dB  bandwidth  and  50  gsec  delay.  Chirp  waveforms  of  very  large 
time-bandwidth  product  (BT<2”“524,288)  were  synthesized 
digitally.  The  delay  range  was  limited  to  36  X  27  /isec  (4.-3  aspect 
ratio)  to  increase  the  signal  bandwidth.  The  image  plane  sampling 
was,  therefore,  380  +  36  /isec  s  10.6  MHz  in  one  delay  dimension 
and  488  +  27  <isec  a  18.1  MHz  in  the  other  dimension.  The  overall 
system  frequency  response  (MTF)  was  approximately  3  dB  lower  at 
the  Nyquist  limit  (2  samples  per  cycle)  5  and  9  MHz,  respectively. 
An  example  ambiguity  function  of  a  short  psuedo-random  code  is 
shown  in  Figure  8.  The  Doppler  range  was  1-3.5  kHz  determined  by 
the  chirp  rate.  The  code  repetition  period  was  6.2  nsec  (31  length,  5 
MHz  rate),  therefore.  4  correlation  peaks  are  evident  in  the  am¬ 
biguity  function;  the  Doppler  was  2  kHz.  The  time  bandwidth  pro¬ 
duct  was  5  MHz  x  33.34  msec  **  166,700.  No  post-detection  pro¬ 
cessing  has  been  applied  to  the  output  video.  A  signal  dependent 
bias  variation,  proportional  to  the  pairwise  cross-correlations  be¬ 
tween  inputs,  has  not  been  filtered.  The  image  has  several 
blemishes  due  to  the  CCD  camera.  The  sensor  dynamic  range  is 
approximately  60  dB.  and  the  maximum  signal-to-bias  ratio  was 
25%:  the  resultant  input  signal  dynamic  range  was,  therefore,  48 
dB 

VI.  SPECTRAL  ANALYSIS 

A  noncoherent  optical  time  integrating  approach  to  two- 
dimensional  complex  spectral  analysis  is  described.  This  method  h 
an  extension  of  the  technique  described  for  one-dimensional 


analysis.  Very  large  time-bandwidth  (greater  than  10s)  spectral 
analysis  of  electrical  signals  can  be  performed  using  this  time 
domain  approach,  without  requiring  storage  of  the  signal.  Spectral 
analysis  of  raster  scanned  video  imagery  can  be  performed  as  well. 

The  optical  system  (Figure  7)  described  in  the  last  sections  is 
used  in  this  example  as  well.  The  algorithm  is  implemented  with 
three  reference  chirp  waveforms.  The  source  is  modulated  with  the 
input  signal,  which  is  premultiplied  by  a  chirp,  and  the  other  chirps 
modulate  A01  and  A02.  The  chirp  rates  are  chosen  such  that  the 
product  of  the  three  chirps  result  in  a  temporal  sine  wave,  with  fre¬ 
quency  that  is  varying  from  pixel  to  pixel  on  the  detector  array. 
Each  detector  element  integrates  the  product  of  this  sine  wave  times 
the  input  signal,  and,  in  this  fashion,  produces  the  spectrum.  The 
frequency  difference  between  output  lines  in  one  dimension  is 
taken  to  be  N  times  the  frequency  difference  between  lines  in  the 
other  dimension,  where  N  is  the  number  of  detectors  per  line.  For  a 
square  array,  the  number  of  spectral  samples  is  proportional  to  N2. 
The  technology  is  currently  limited  by  detector  array  size,  and 
detector  readout  rate. 

A  two-dimensional  chirp  algorithm  with  three  complex 
reference  chirps  is  described  next.  The  real  implementation  is  per¬ 
formed  with  frequency  offset  chirps,  similar  to  the  previous 
examples.  Define  the  transform  (optical  processor  output)  by 

iNT/2 

s(t)a(t)b(t  —  r,  )c(t  —  r2)dt  (20) 

-NT/2 


where 

alt)  =  2^e')a*(,‘nT)1  rect 
n 

b(t)  -  e  )°»,vN  rect(t/NT) 

c*(t)  -  a(t)b(t)  .  (21) 

Waveforms  aft)  and  eft)  are  periodic  chirps  with  period  T  and  chirp 
rates  a  and  a(N  — 1)/N,  respectively.  Waveform  b(t)  is  a  low  rate 
chirp  with  period  NT  and  rate  a/N.  The  chirp  rates  are  chosen  such 
that  a  constant  difference  frequency,  f  —  a[r,  +  (N  —  l)rj|/N,  is 
generated  by  the  product  a(t)b(t  —  r,  )c(t  —  t2).  Frequency  discon¬ 
tinuities  occur  during  intervals  nT  <  t  <  nT  +  t,  due  to  the  chirp 
transition  traversing  the  acoustic  delay  aperture.  This  effectively 
reduces  the  integration  time  by  the  factor  (T-  r,)/T.  However,  this 
may  be  completely  compensated  by  increasing  the  chirp  band¬ 
width.  This  effect  is  ignored  in  the  following  analysis,  thus 
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a(t)b(t~Ti)c(t-r,)  • 

jUwomTri  -2wat(T,  +  (N-l)r,|/N  +  $<r,.T,)| 
n 

>rect(ifL)rec,(^)  (22> 

where 

0<r„r j)  =  or[r{  +  (N-l)rj]  /n  . 

The  spectral  phase  weighting,  may  be  compensated  post¬ 

detection.  Equation  (20)  becomes 

R(t,,t2)  =  e^^NT  sinc(Tf)  Sff-afa  +  fN-UrJ/N] 

•  ^sinc[NT(f  +  arj  — n/T]df  .  (23) 

n 

As  evidenced  by  the  term  S(a(ri  +  (N  —  1)ti)/N1,  the  variables  r, 
and  T]  may  be  interpreted  as  fine  and  coarse  frequency,  respec¬ 
tively.  Integration  over  N  chirp  repetitions  has  created  a  comb  sam¬ 
pling  along  the  coarse  frequency  axis,  with  comb  samples  a  func¬ 
tion  of  the  fine  frequency  axis.  The  spectral  resolution  is  propor¬ 
tional  to  1/NT  where  NT  is  the  total  integration  period. 

VII.  SUMMARY 

Time  integrating  optical  techniques  for  one-  and  two-dimensional 
complex  signal  processing  have  been  described.  Signal  processing 
architectures  utilizing  actively  generated  reference  waveforms  such 
as  chirps  realize  a  wide  range  of  algorithms  and  variable  time  inte¬ 
gration  allows  further  flexibility.  Particularly  attractive  is  the 
multipurpose  capability  of  the  time  integrating  optical  processor, 
e.g.,  the  same  optical  system  is  used  for  both  spectral  analysis  and 
ambiguity  function  processing. 

The  technique  of  interferometric  detection  with  electronic 
reference  was  described.  In  this  implementation,  a  local  oscillator  is 
added  to  the  acoustic  modulation  rather  than  a  coherent  optical 
reference  added  to  the  image.  This  method,  therefore,  permits  non¬ 
coherent  optical  implementation,  with  the  advantages  of  directly 
modulated  diode  light  sources  and  increased  immunity  to  artifacts 
of  coherent  optical  imaging  systems.  The  correlation  is  performed 
at  a  carrier,  thereby  enabling  complex  operation.  Furthermore,  the 
interferometric  method  circumvents  the  difficult  requirement  for 
acoustic  modulation  at  a  high  diffraction  efficiency  bias  point,  that 
is  necessary  for  linear  operation  in  the  noninterferometric  techni¬ 
que 

The  optical  implementation,  utilizing  acousto-optic  input  and 
integrating  image  sensor  output  devices,  creates  a  processor  that  is 
highly  compatible  with  signal  processing  systems.  Such  implemen¬ 
tation  affords  very  large  time-bandwidth  product  signal  processing 
without  the  input  signal  storage  requirement  associated  with  spatial 
integrating  methods.  The  requirement  for  high  resolution,  large 
dynamic  range  output  image  sensors  becomes  the  key  device  limita¬ 
tion.  Acousto-optic  devices  are  available  with  time-bandwidth  pro¬ 
duct  much  greater  than  the  number  of  resolvable  image  samples. 


The  key  attributes  of  time  integrating  techniques  are  summa¬ 
rized  in  the  following:  extremely  large  time-bandwidth  correlations 
may  be  performed,  independent  of  the  device  time-bandwidth  pro¬ 
duct:  a  flexible,  multipurpose  processor  is  realized  by  use  of  ac¬ 
tively  generated  reference  waveforms  and  variable  time  integra¬ 
tion;  an  important  class  of  two-dimensional  algorithms,  including 
complex  spectral  analysis  and  ambiguity  function  processing,  may 
be  performed  without  having  to  store  the  entire  time  history  as  a 
spatial  record;  the  system  may  be  implemented  with  noncoherent 
diode  sources,  acousto-optic  devices,  and  integrating  image 
sensors. 
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Abstract 


Utilization  of  both  space  and  time  integration  in  acousto  optic  signal  processing  over¬ 
comes  some  of  the  limitations  of  purely  space  integrating  and  purely  time  integrating 
techniques.  A  one-dimensional  configuration  for  coherent  hybrid  spectrum  analysis  is  used 
to  show  that  this  approach  exhibits  superior  dynamic  range  characteristics,  including  a 
unique  sidelobe  suppression  property.  A  coherent  hybrid  ambiguity  function  processor  is 
discussed  to  illustrate  the  important  advantages  that  may  be  exploited  by  using  both 
space  and  time  integration  in  delay-Doppler  plane  processing. 

Introduction 

Acousto-optic  techniques  are  becoming  increasingly  attractive  for  wide  band  parallel 
processing  because  of  the  real-time  capability,  the  maturing  system  architectures  and 
components,  and  the  potential  for  small  size  and  low  power  requirements.  Time  integration 
techniques  1 “ 4  in  one  and  two  dimensional  systems  have  overcome  some  of  the  device-related 
limitations  by  permitting  long  time  processing  for  spectrum  analysis,  correlation, 
ambiguity  function  processing  and  other  two-parameter  processes.  Hybrid  space/time  inte¬ 
gration  approaches  6  have  been  shown  to  retain  some  important  characteristics  of  space 
integration  and  permit  extended  time  processing  by  time  integration.  This  paper  presents 
a  summary  of  some  recent  results  on  hybrid  space/time  integration  approaches  to  signal 
processing.  Two  topics  are  addressed.  First,  some  further  considerations  of  coherent 
hybrid  spectrum  analysis  are  presented,  emphasizing  aspects  of  system  dynamic  range.  It 
is  shown  that  the  combination  of  space  and  time  integration  leads  to  a  superior  performance 
that  makes  this  approach  attractive  even  for  one-dimensionc 1  systems.  Most  noteworthy  is 
the  unique  behavior  of  hybrid  systems  of  suppressing  sidelobes. 

The  second  part  of  this  paper  extends  the  coherent  hybrid  approach  to  the  generation 
of  cross  ambiguity  functions.  It  is  shown  that  the  combination  of  space  and  time  integra¬ 
tion  here  also  leads  to  superior  performance  in  dynamic  range  and  signal  handling  ability. 

Spectrum  Analysis 

A  relatively  simple  but  highly  successful  application  of  acousto  optics  in  signal  pro.- 
cessing  is  the  familiar  acousto  optic  power  spectrum  analyzer.  This  is  referred  to  as 
a  space  integrating  technique  because  the  Fourier  transform  integral  is  over  a  spatial 
coordinate.  The  usefulness  of  such  a  technique  derives  from  the  fact  that  processing  of 
many  channels  is  performed  in  parallel  over  a  wide  bandwidth.  Important  limitations 
encountered  with  this  approach  are  that  the  frequency  resolution  is  limited  to  the  inverse 
of  the  acoustic  transit  time,  the  number  of  channels  is  limited  to  the  time-bandwidth  pro¬ 
duct  of  the  device,  and  the  dynamic  range  is  frequently  severely  limited  by  the  detector 
noise,  since  the  detector  senses  the  signal  power. 

Time-integrating  AO  techniques  overcome  these  limitations  by  coherent  detection  of  the 
signal  by  many  local  oscillators  distributed  over  the  signal  band.  Here  the  complex 
spectral  amplitude  is  detected  so  that  the  detector  noise  is  significantly  less  important. 
The  frequency  resolution  is  the  inverse  of  the  integration  time  and  is  (along  with  the 
analysis  bandwidth)  electronically  scaleable.  The  number  of  parallel  channels  can  be 
extended  in  two  dimensions  to  the  product  of  the  time-bandwidth  products  of  two  Bragg 
cells.  Time-integrating  techniques,  however,  have  their  own  limitations,  including  the 
constraint  to  lower  bandwidths  and  a  dynamic  range  that  degrades  as  the  number  of 
occupied  channels  increases. 

Previous  studies  '  have  shown  that  hybrid  configurations  utilizing  both  space  and 
time  integration  can  retain  important  performance  characteristics  of  both  space  and  time 
integrating  approaches.  Here  we  continue  the  discussion  of  hybrid  spectrum  analysis  by 
describing  a  one  dimensional  configuration,  focusing  on  the  characteristics  that  princi¬ 
pally  affect  the  system's  dynamic  range.  Such  a  system  is  shown  schematically  in  Figure  1. 
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Figure  1.  One  dimensional  coherent  hybrid  spectrum  analyzer 

Coherent  light  is  diffracted  from  two  Bragg  cells,  one  driven  by  the  signal  and  the 
other  by  a  periodic  chirp,  such  that  they  are  both  Doppler-shi f ted  in  the  same  direction. 
Both  optical  waves  are  Fourier  transformed  by  a  lens  onto  a  detector  array,  where  they 
interfere  to  form  a  fringe  complex.  This  fringe  pattern  is  integrated  over  time  T.  equal 
to  the  Bragg  cell  delay  time.  The  resulting  exposure  is  a  sampling  of  the  signal's  am- 
ilitude  spectrum  on  a  spatial  carrier  determined  by  the  mean  relative  angle  between  the 
incident  waves. 

The  signal  in  Figure  1  is  assumed  to  lie  in  the  band  of  the  Bragg  cell.  If  v  is  the 
acoustic  velocity,  x  is  the  coordinate  in  the  detector  plane  and  z  is  the  coordinate 
ir  the  Bragg  cell  driven  by  signal  s  t) ,  the  amplitude  at  the  detector  is 


S(x,t) 


z/v) 


-  i  2 1  ax  z 
e 


dz 


f  i2xf*t 
Je  u 


(f ') 


sine  TTT^avx-f')  df' 


(1) 


where  u(f)  is  the  signal  amplitude  spectrum,  a  =  1/XF,  X  is  the  optical  wavelength  and  F 
is  the  focal  length  of  the  lens.  Consider  a  CW  signal  of  frequency  f.  Then 

S(x,t)  =  SQe  sine  nT^lavx-f).  (2) 


The  periodic  chirp,  unlike  the  time-integrating  techniques,  has  a  chirp  length  and 
period  equal  to  the  Bragg  cell  delay  T..  Its  bandwidth  represents  the  system’s  analysis 
bandwidth  and  can  be  as  large  as  that  of  the  Bragg  cell.  The  Fourier  transform  of  the 
light  diffracted  from  this  periodic  chirp  can  be  considered  a  local  oscillator  in  the 
detector  plane  that  has  a  frequency  that  is  a  linear  function  of  position  on  the  array. 

It  can  be  approximately  described,  as  shown  in  the  Appendix  by 

R(x,t)  =  Ro  £  ei2iipt/T1  +  in<p>sinc  tiTj  [avx  -  p/TJ  (3) 

The  normalized  exposure  after  some  time  T2  >  T^  is 


E(x)  =  i 


/ 


T2/2 


2  -T2/2 


lRol  + 


|R(x,t)elkx  +  S (x, t) | 2  dt 


2  2 

|SQ|  sine  nTjLavx  -  f) 


2,^  sinctrTjLf  -  p/T^  sincsTj^Lavx  -  p/Tj) 


x  sincnT^  (avx  -  f)  cos(kx  +  ip  +  Slip)). 


(4) 


SPIE  Vol  232  1980  Intarnational  Optica!  Computing  Confaranca ! 1980 }  83 


I 


! 


I’  t 


When  Tj  =  it  can  be  shown  that  the  Whittaker-Shannon  sampling  theorem  leads  to  the 
identity 


sincuT^ [f-p/T^l  sincnT^ [avx-p/T^ ] 


=  ei«(avxT1)+i21-voavx  sinc7lTi  ,f _avx,  (5) 


where  vc  =  f  and  Af  are  the  minimum  frequency  and  width  of  the  band. 

Substitution  of  Equation  S'  into  Equation  4  yields  the  signal  term 

E  (x)  =  2|R  S  |  sinc^irT,  (avx  -  f)  cos(kx  +  ij,  +  f}*(ovxT, )  ]  .  (6) 

S  O  O  1  1 

2 

Thus  the  signal  envelope,  the  peak  of  which  is  proportional  to  |S  |,  has  a  sine  form, 
rather  than  sine  as  in  the  time-integrating  and  space  integrating°techniques .  An  important 
consequence  of  this  is  that  the  sidelobes  are  down  by  a  factor  of  two  in  dB  compared  with 
those  techniques.  Sidelobes  of  the  very  bright  zero  order  (undiffracted  by  the  acoustic 
wave)  that  extend  into  the  signal  band  are  similarly  suppressed.  This  inherent  sidelobe 
suppression  is  the  result  of  the  combined  effect  of  the  localization  of  the  signal  illumi¬ 
nation  at  the  detector  by  spatial  integration  (i.e.,  the  optical  Fourier  transform)  and 
localization  of  the  interference  fringe  pattern  by  time  integration.  In  real  space 
integrating  and  time  integrating  systems  the  profile  of  a  spectral  line  is  not  actually  a 
sine  function  but  is  modified  by  the  apodization  of  the  optical  wave  illuminating  the 
Bragg  cell,  among  other  things.  Similar  apodizing  behavior  applies  to  hybrid  systems; 
moreover,  if  P(avx-f)  is  the  profile  of  the  spectral  amplitude  for  a  spape  or  time  inte¬ 
grating  spectrum  analyzer,  that  for  a  hybrid  system  can  be  shown  to  be  P/(avx-f). 

A  source  of  noise  in  optical  processors  that  can  limit  the  dynamic  range  is  light 
scattered  from  optical  elements,  particularly  the  Bragg  cells.  Such  scattered  light  is 
not  Doppler  shifted,  and  therefore  in  time  integrating  and  hybrid  systems,  the  only  terms 

in  the  exposure  that  do  not  integrate  to  zero  are  products  containing  only  scatter 

elements.  If  sn(x)  and  r  (x) exp ( i**)  are  the  scatter  amplitude  components  at  the  detector, 

the  exposure  term  2 1  sn  r^Jcos  (kx+41)  has  a  smaller  variance  than  the  signal  power  by  a 

"”5  2 

factor  |s  rn|  /  |R  S|  .  Therefore  the  dynamic  range  referred  to  scatter  noise  is 
twice  as  large  (in^dB)  for  systems  using  time  integration  than  for  purely  space  integrating 


systems.  sue  i  °n  j  i  - i  ~ —  *  ■»—  — — —  ~r - 

frequencies,  outnof  the  signal  band) . 

It  will  be  noted  that  in  Eq.  4  the  bias  term  contributed  by  the  signal  is  localized 
and  does  not  contribute  a  bias  to  other  signals  in  the  band.  Unlike  time  integrating 
techniques  that  exhibit  reduced  dynamic  range  for  multiple  signals  because  of  a  uniform 
bias  contributed  by  each  signal,  hybrid  techniques  retain  a  high  multiple  signal  or  broad 
band  capability. 

Let  us  now  consider  the  dynamic  range  as  limited  by  detector  noise.  The  output  voltage 
of  a  linear  detector  array  is  proportional  to  the  exposure,  where  we  can  write 


terms  are  also  small,  but  are  furthermore  at  low  spatial 


V(x)  = 


2a|RS|  cosi)i(x) 


The  signal  term  remaining  after  bandpass  filtering  is 

Vg(x)  =  2ajRS|  cosiHx)  . 

?  2 

This  is  a  maximum  when  at | R |  =  u|S  |  =  v_  t/4 ,  where  V  .  is  the  saturation  voltage. 

Thus  the  maximum  detected  signal  power  is  a  8  r 


1 2aRS^cos$ (x) 


define  the  minimum  detected  signal  power  to  be 


V  sat/8  • 


| 2aRS  cost ( x )  I  2  =  o2 

m  • 

vhere  0  is  the  rms  detector  noise  over  one  frequency  resolution  element.  The  dynamic 
range  related  to  detector  noise  is  therefore 


IVSml 


(Vsat/'’>  /8 


84  /  SPIE  Vo!  232  1980  Intarnational  Optical  Computing  Conlaranca  11980) 


COHERENT  HYBRIO  OPTICAL  PROCESSORS 


A  space  integrating  power  spectrum  analyzer  has  a  dynamic  range  DR  :  V  t/a  30  that 
the  advantage  in  dynamic  range  of  interferometric  techniques  is  described  By  the  relation- 

ShlF  DR  =  (DRq)2/8  . 

This  advantage  is  directly  analogous  to  the  well-known  improvement  in  dynamic  range  of 
optical  heterodyne  detection  compared  with  direct  detection.' 

Ambiguity  Function  Processing 

Two  dimensional  AO  techniques  have  been  found  useful  for  simultaneous  display  of  time 
delay  and  Doppler  shift  of  two  signals  in  the  form  of  their  cross-ambiguity  function 


Xlx.f!  = 


=/Sl<t>S2 


dt 


There  are  numerous  implementations  2~ “  of  time  integrating  techniques,  both  coherent  and 
incoherent,  all  of  which  are  characteri zed  by  the  property  that  the  detector  is  in  an 
image  plane  of  one  or  both  of  the  Bragg  cells  driven  by  the  signals.  Let  T.  be  the 
Bragg  cell  delay.  The  exposure  of  the  two  dimensional  detector  with  coordinates  t  and  f 
after  time  T  is  generally  of  the  form 


where  signifies  convolution,  B  is  a  bias  term  that  is  more  or  less  uniform  over  the 
detector,  C  is  a  zero  mean  correlation  noise,  the  rect  function  reflects  the  delay  range 
limitation  to  the  Bragg  cell  delay  (c  =  1  or  2,  depending  on  the  particular  implementation) 
and  the  sine  function  indicates  that  the  Doppler  Rayleigh  resolution  is  the  inverse  of 
the  integration  time.  The  cos  term  represents  a  spatial  carrier  that  has  been  arbitrarily 
chosen  to  be  in  the  i  direction. 


Let  us  consider  as  an  example  a  long  pseudorandom  signal  of  bandwidth  Afs.  The  cross 
ambiguity  function  between  the  signal  and  a  shifted  and  delayed  version  over  a  time  window 
T,  is  shown  schematically  in  Figure  2,  exhibiting’  the  familiar  "thumbtack” 
shape  with  a  sharp  peak  of  width  1/Afs  and  I/T2  in  the  delay  and  Doppler  direction  respec¬ 
tively,  (when  the  peak  is  not  too  far  from  the  center).  The  correlation  noise  or 
self-clutter  is  spread  over  a  delay  and  Doppler  area  &fg  T2.  Figure  3  shows  the  cross¬ 
ambiguity  function  measured  by  a  time  integrating  processor  integrating  for  time  T,.  It 
is  a  small  area  cut  out  of  the  function  of  Fig.  2,  where  the  delay  range  is  limited  by  the 
Bragg  cell  delay,  according  to  Eq.  7,  and  the  Doppler  range  is  chosen  for  convenience. 

The  bias  term  is  also  included.  A  point  to  be  noted  is  that  if  one  of  the  two  signals  to 
be  processed  contains  more  than  one  delayed  or  shifted  component,  each  component  adds  a 
bias  contribution,  thus  reducing  the  dynamic  range  per  component  in  a  manner  analogous  to 
the  multiple  signal  dynamic  range  reduction  in  time  integrating  spectrum  analyzers. 


Figure  2.  Cross  Ambiguity  Function  of  a  Pseudorandom  Signal 
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Figure  3.  Cross  Ambiguity  Function  for 
Time  Integrating  Processor 

Space  integrating  techniques  *  have  been  used  to  generate  the  cross  ambiguity  func¬ 
tion;  the  constraint  encountered  when  acousto-optic  devices  are  used  in  real  time  is  that 
the  Doppler  resolution  is  limited  to  the  inverse  of  the  Bragg  cell  delay  time.  Further¬ 
more,  unless  the  optical  wave  is  coherently  detected  the  technique  suffers  from  the  same 
dynamic  range  limit  as  a  power  spectrum  analyzer;  that  is,  the  absolute  square  of  the 
ambiguity  function  is  detected,  permitting  the  detector  noise  to  be  a  more  serious 
problem. 

The  combination  of  space  and  time  integration  can  be  applied  to  ambiguity  function 
processing  to  provide  improved  performance.  Consider  the  optical  schematic  of  Figure  4. 
It  consists  basically  of  a  space  integrating  ambiguity  function  generator  designed  to 
facilitate  expansion  of  the  Doppler  coordinate  so  that  a  small  portion  of  the  Doppler 
range  fills  the  detector  array  and  a  spatially  distributed  local  oscillator  that  provides 
a  coherent  reference  for  time  integration.  The  first  Bragg  cell  driven  by  Sj(t)  is 
imaged  with  magnification  m.  onto  the  second  cell  orthogonal  to  it,  driven  by  S2<t). 
Diffraction  orders  of  the  two  cells  of  opposite  sign  are  selected.  This  image  plane  is 
rotated  by  angle  0  with  respect  to  coordinates  x  and  y  with  respect  to  which  Iburier 
transformation  and  imaging  with  magnification  m2,  respectively,  are  performed.  The 
resulting  optical  amplitude  at  the  detector  is 

A(t 


,  f ;  t )  =  ce 


-i<l/(T,f) 


rect  r 
2T, 


j*(T,f)  e  i2,,f(t-T/2)|*P(T,f)  , 


(8) 


i 


) 


where 

t  =  y  cosa  /vm^m^  ,  f  =  xv/(XF  cos3) 

are  the  delay  and  Doppler  coordinates,  respectively,  F  is  the  focal  length  of  the  trans¬ 
form  lens,  and  where 

P  ( T ,  f )  =  (T1-  |  t  | )  sine*  (T  -  |  t  | )  f 

represents  the  Doppler  resolution  function.  The  best  resolution  is  seen  to  be  1/T^  cor¬ 
responding  to  small  delays. 

The  objective  of  the  hybrid  scheme  is  to  cut  out  a  small  band  of  the  ambiguity  function, 
reducing  the  Doppler  range,  and  improve  the  Doppler  resolution  by  time  integration. 

Adding  the  reference  wave  R  exp  iC2n (v/AF) xt  -  kx]  and  integrating  the  absolute  square  of 
the  result  over  time  T2  »T^  leads  to  an  exposure 

E(i,f)  =  R2  +  jcITj-M)  x  <  T  ,  f )  |  2  #  sind2nT^f 


for  |  t  |_^  Tj 


2R [Tj- | t | ] |cx ( i , f ) *  sincrT^f I 


cos (k  ’x+i^  ’ ) 


(9) 
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Figure  4.  Schematic  of  Coherent 
Hybrid  Ambiguity  Function  Processor 


For  our  psendorandom  signal  this  is  represented  in  Figure  5.  The  broad  Doppler  peak 
has  been  reduced  by  time  integration  to  a  peak  with  a  spatial  carrier  that  is  1/T2  wide. 
An  important  characteristic  of  this  measured  cross  ambiguity  function  is  that  the  bias 
contribution  from  the  signal  arm  of  the  interferometer  is  confined  to  the  delay  bias. 
Multiple  signal  components  generating  peaks  separated  by  at  least  the  delay  resolution 
1 / A f s  do  not  degrade  the  dynamic  range  as  a  result  of  bias  accumulation,  as  is  the 
case  for  purely  time  integrating  processors.  If  the  chosen  Doppler  range  is  larger  than 
the  Doppler  resolution  of  the  space-integrating  processor  above,  then  some  localization 
of  bias  in  the  Doppler  direction  is  also  maintained. 

Because  of  the  time  integration  aspect,  optical  scatter  noise,  as  in  the  spectrum 
analyzer,  integrates  out  except  for  small  second  order  terms. 


Figure  5.  Ambiguity  Function  For  Coherent  Hybrid  Processor. 
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Summary 

Using  the  examples  of  one  dimensional  spectrum  analysis  and  two  dimensional  ambiguity 
function  processing,  some  of  the  limitations  of  purely  space  integrating  and  purely  time 
integrating  techniques  for  acoustooptic  signal  processing  can  be  overcome  by  utilizing 
both  space  and  time  integration.  Dynamic  range  advantages  resulting  from  the  coherent 
detection  aspect  of  interferometric  techniques  are  enhanced  by  optical  scatter  noise 
reduction,  localized  bias  for  superior  multiple  signal  handling  ability,  and,  in  the  case 
of  spectrum  analysis,  a  reduction  in  spectral  sidelobes.  The  inherent  sidelobe  suppression 
property  of  coherent  hybrid  spectrum  analysis  is  an  important  factor  in  many  applications, 
particularly  those  involving  threshold  detection,  and  is  characteristic  of  no  other  known 
optical  technique. 

The  author  acknowledges  support  for  this  work  by  the  Air  Force  Office  of  Scientific 
Research  under  contract  F49620-78-C-0102. 


Appendix 


An  acoustooptic  device  having  sound  velocity  v  is  driven  by  a  periodic  chirp  of  length 
T ^ ,  period  Td,  bandwidth  Af  and  minimum  frequency  fm  such  that  the  amplitude  of  the 
up-shifted  diffracted  light  is 


r  ( t ,  f, )  =  J  eila(t-nT1+C/v)2  +  b(t-nT1+C/v)]  rect  <t-nT1+C/v) /Tl  recU/D 


(A- 1) 


where 


it  Af 


b  =  2nf 


At  the  detector  plane,  an  optical  Fourier  transform  plane  of  the  Bragg  cell,  the 
convolution  theorem  yields 


R(x,t) 


,t)  sincnaDx’  dx' 


where 


/*  OD  ja  OO 

r(t,£)  exp(-i2na£.x)  dC  =  D  /  V(x-x‘, 

—  oo  —  ao 

_  oo 

v(x,t)  =  2  /e1'a(t-nTl+^v)  +  b(t-nT1+C/v)  ]  -  l2l’axtrect(t_nTi+C/v)/Ti  dr  (A. 


( A-2 ) 


3) 


The  change  of  variables  z  =  t-nT-^+C/v  leads  to 


OO 

V(x.t)  =  v  2  .-i2«avx(t-nT1)  f 
n  =  -a>  J-c 


Using  the  identity 
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e  rect(z/T.)  dz  . 
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(A- 5) 


where  6(x)  is  the  Dirac  delta  function,  we  have 


V  (x,  t)  =  c  e 


i2favxt 


6  lx-  (p/avT 


oo 

1>> 


Iaz  -  (2navx-b) z) 


rectlz/Tj^)  dz 


( A-6 ) 


Substitution  of  Equation  A-6  into  Equation  A-2  yields 


R(x, t) 


C1  ^*(P/‘T1>  el2T,pt/Tl 
P 


sine  nT^ [ avx-p/T^ ] 


(A-7) 
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where  is  a  constant  and 


■Mp/T1) 


(b-2np/T1) z) 


rect ( z/T^) 


dz 


(A-8) 


The  function  1 4> ( p/T -) J  describes  the  uniformity  of  the  modulus  in  the  transform  plane 
of  the  Bragg  cell  in^our  model.  For  convenience  of  discussion  of  spatially  distributed 
local  osollators  we  can  derive  an  approximation  to  Equation  A-8  in  the  limit  of  large 
N.  (/sing  the  convolution  theorem, 


where 


GKp/T^-C) 


sine  tiT^C1  d£’ 


(A- 9) 


GNp/V-C’J  =  c2  e‘i"2l(P/T1)‘(b/Zl,>-S']2/a  .  (A_10 

2 

In  the  expansion  of  the  exponent  in  Equation  A-10  the  (£')  term  can  be  considered 
negligible,  because  at  the  first  null  of  the  sine  function  in  Equation  A-9,  £'=1/^  and 


■n2(£')2/a  =  "/AfT,  =  n/N  «1  . 


Therefore  we  can  write  approximately 


R(x,  tj 


Jei21,pt/T1 

P 


l  [  (2np/Tj )  -b]  /4a  sinc  „Tijqvx  _ 

(A- 11) 


where  p  runs  over  intergers  from  f  I.  to  (f  +Af)T. ,  and  where  we  have  observed  that 
D  =  vT^.  Cj  is  a  complex  constant.  m 
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TIME-INTEGRATING  PROCESSORS  USING  BRAGG  CELLS 


Todd  R.  Bader 
ESL  Incorporated 
Sunnyvale,  CA 


Acousto-optic  methods  are  currently  gaining  wide  use  in 
electronic  signal  processing  systems  in  a  variety  of  appli¬ 
cations,  particularly  communications  and  radar  signal 
processing.  The  motivation  for  acousto-optic  techniques  is  the 
associated  wide  bandwidth  and  parallel  processing  capability. 

The  blossoming  realization  of  high  performance  hardware  is  a 
result  of  the  currently  maturing  technologies  of  key  devices, 
including  the  acousto-optic  devices  themselves  and,  equally 
important,  coherent  and  noncoherent  light  sources  and  high 
quality  optical  detectors. 

The  acousto-optic  devices  (Bragg  cells)  used  in  signal 
processing  systems  are  real-time  spatial  light  modulators;  that 
is,  a  signal  impressed  on  the  piezoelectric  transducer  is 
represented  (due  to  an  acousto-optic  interaction  in  the  medium) 
as  a  traveling  spatial  variation  of  the  light  amplitude 
emerging  from  the  device.  The  light  diffracted  by  a  component 
of  this  traveling  sound  wave  is  also  Doppler  shifted  in 
optical  frequency  by  the  frequency  of  that  component,  a  signif¬ 
icant  property  relative  to  the  use  of  A-0  devices  for  time- 
integrating  coherent  processing. 

Bragg  cells  are  currently  available  in  high  quality 
devices  operating  in  bandwidths  that  range  from  several  MHz  to 
perhaps  2  GHz,  and  with  transit  time  X  bandwidth  products  up 
to  about  2000.  The  required  bandwidth  determines  the  acousto¬ 
optic  material  used.  Generally  the  materials  with  the  best 
acousto-optic  coupling  have  the  lowest  bandwidth  capability. 

The  most  straightforward  application  of  acousto-optic 
devices  for  signal  processing  is  a  power  spectrum  analyzer  (or 
Bragg  receiver)  in  which  an  array  of  photodetectors  samples 
the  optical  intensity  pattern  in  the  back  focal’  plane  of  a  lens 
that  Fourier  Transforms  the  coherent  light  amplitude  diffracted 
from  a  Bragg  cell.  AIbo,  a  variety  of  simple  optical  systems, 
both  coherent  and  noncoherent,  can  be  implemented  to  measure 
cross-correlation  functions  between  two  signals.  Bragg 
receivers  and  time-integrating  optical  correlators  are  partic¬ 
ularly  useful  for  sensitive  detection  of  narrowband  and  spread 
spectrum  signals,  respectively.  Both  algorithms  use  a  time- 
integrated  detected  signal  to  improve  sensitivity.  This 
post-detection  integration  has  the  additional  benefit  of 
reducing  the  processor  output  data  rate,  easing  the  requirements 
placed  on  the  detector  array.  Interferometric  techniques  for 
coherent  detection  may  be  used  for  complex  amplitude  spectrum 
analysis  at  the  cost  of  more  sophisticated  detector  array  and 
post-detection  processing  technology. 
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Two-dimensional  systems  composed  of  a  light  source,  Bragg 
cells,  two-dimensional  detector  arrays  and  simple  passive 
components  can  be  configured  to  perform  a  variety  of  processing 
algorithms.  High  resolution  spectrum  analysis  is  performed  by 
extending  frequency  resolution  in  the  transverse  direction  of 
a  coherent  spectrum  analyzer  by  time-integration.  Complex 
spectra  of  imagery,  for  example,  can  be  generated  in  this  way 
by  providing  an  input  electronic  signal  in  the  form  of  a  raster 
scanned  image.  Real-time  cross  correlation  of  two  images  is 
performed  by  generating  an  inverse  2-D  transform  of  the  product 
of  their  spectra.  Another  2-D  process  well-matched  to  A-0 
techniques  is  the  generation  of  ambiguity  functions.  This  is  a 
valuable  method  of  measuring  the  cross  correlation  of  a  signal 
with  a  version  of  it  having  an  unknown  Doppler-shif t . 

Acousto-optic  processing  performance  is  presently  limited 
by  detector  array  technology.  Although  systems  are  currently 
being  fielded  with  high  performance  A-0  processors,  these  are 
mostly  one-dimensional.  Detector  technology  is  rapidly 
advancing,  however,  and  the  near  future  should  see  a  growing 
volume  and  variety  of  one-  and  two-dimensional  acousto-optic 
signal  processors  for  numerous  applications. 
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